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1.  BACKGROUND 


In  most  instances,  Army  Regulation  (AR)  385-64  requires  that  quantity-distance  computations  to 
determine  ammunition  storage  siting  limitations  be  based  on  the  total  weight  of  all  mass-detonating 
explosives  at  the  site.  This  requirement  is  made  under  the  assumption  that  the  blast  waves  emanating  from 
storage  subdivisions  (ammunition  "stacks")  as  they  sequentially  detonate  will  coalesce  with  one  another 
before  their  amplitudes  become  negligible.  If  it  can  be  shown  that  barriers  between  stacks  prevent  or 
sufficiently  delay  propagation  of  detonation  between  them,  distances  may  be  determined  using  the 
explosive  weight  of  each  stack  individually.  The  regulation  provides  a  crude  criterion,  based  on  initiation 
delay,  for  determining  when  coalescence  will  occur  for  the  explosion  of  two  closely  spaced  stacks.  This 
criterion  has  not  been  verified  for  complicated  ammunition  storage  arrangements  and  may  be  too 
restrictive,  particularly  where  significant  spatial  separation  of  stacks  exists. 

Experience  with  the  quantity-distance  requirements  has  shown  that  they  are  often  difficult  to  meet  and 
waivers  of  the  regulations  have  been  frequently  requested.  If  the  actual  coalescence  pattern  from  a 
complex  array  of  ammunition  stacks  could  be  estimated,  a  more  accurate  assessment  of  the  hazards 
produced  could  be  made,  and  significant  reductions  in  real-estate  requirements  might  result  In  principle, 
this  could  be  accomplished  by  numerical  simulation  of  the  blast  environment  in  the  region  of  interest. 
However,  tWs  approach  requires  immense  computer  resources,  even  for  the  simplest  arrangements,  and 
is  beyond  the  reach  of  the  typical  user.  Alternatively,  approaches  based  on  simplified  analyses  which  have 
been  developed  in  the  nuclear  blast  community  could  be  applied.  Specifically,  algorithms  found  in  the 
DNA  Nuclear  Blast  Standard  (1  KT)  (Needham  and  Crepeau  1981),  a  computer  program  which  describes 
the  blast  environment  produced  by  a  1-KT  nuclear  explosion,  and  the  Low  Altitude  Multiple  Burst 
(LAMB)  model,  a  computer  program  which  describes  the  blast  environment  produced  by  multiple  nuclear 
explosions  of  arbitrary  yield,  served  as  the  starting  point  for  the  present  approach. 

2.  APPROACH 

We  originally  intended  to  adapt  these  programs  to  our  present  purposes.  However,  direct  use  of 
LAMB  required  the  performance  of  many  unnecessary,  time-consuming  computations.  Therefore,  we 
developed  a  new  application  called  BWACO,  intended  to  estimate  pertinent  aspects  of  the  blast 
environment  associated  with  sequentially  detonating,  spatially  distributed  ammunition  stacks. 
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In  its  initial  implementation,  BWACO  made  significant  use  of  1-KT  Standard  and  LAMB  algorithms. 
However,  as  we  investigated  the  performance  of  the  code,  these  were  gradually  abandoned  in  favor  of  new 
algorithms.  The  main  development  effort  was  conducted  using  ARL’s  Cray  X/MP  48  because  graphics 
software  was  readily  available.  However,  once  a  FORTRAN-callable  graphics  package  was  obtained,  the 
code  was  converted  for  use  on  IBM  Personal  Computers  (PCs)  or  compatible  machines. 

This  report  explains  the  assumptions  used  and  documents  the  evolution  of  BWACO  on  the  Cray 
following  its  initial  implementation.  Features  added  in  the  PC  version  are  illustrated  through  application 
to  a  number  of  problems  representative  of  typical  ammunition  storage  configurations.  A  detailed 
description  of  the  code  (including  subroutine  listings)  and  the  initial  versions  of  the  BWACO  User’s 
Manual  are  presented  in  the  appendices. 

3.  BWACO  ASSUMPTIONS 

Sequential  detonations  are  assumed  to  occur  as  a  result  of  propagation  of  detonation  from  one  stack 
to  another.  The  first  stack  of  any  such  pair  to  detonate  is  called  the  donor  and  the  second  is  called  the 
acceptor.  The  first  of  all  the  stacks  to  detonate  is  referred  to  as  the  initial  donor.  Detonation  propagation 
is  assumed  to  occur  at  a  constant  velocity  which  is  a  property  of  the  donor.  This  is  consistent  with  a 
fragment  impact  mechanism  at  least  over  distances  short  enough  to  preclude  significant  deceleration  of 
the  fragments.  The  velocity  provided  must  be  that  of  the  fastest  fragment  produced  by  munitions  in  the 
stack  augmented  to  account  for  focussing  effects.  The  fact  that  no  consideration  has  been  given  to 
aerodynamic  fragment  decelerations  or  to  additional  delays  associated  with  the  penetration  of  barriers 
makes  this  a  worst-case  assumption  (when  sufficiently  high  fragment  velocities  have  been  specified). 

Classical  blast  wave  structure  is  assumed.  The  blast  wave  consists  of  a  shock  wave  and  positive 
overpressure  phase  followed  by  a  negative  overpressure  phase  (when  sufficiently  far  from  the  source). 
The  usual  scaling  laws  are  assumed  to  apply  and  are  used  to  determine  required  blast  wave  characteristics 
with  reference  to  a  standard.  BWACO  requires  shock  and  zero-overpressure  arrival  times  at  specified 
stations.  The  1-KT  Standard  gives  the  shock  and  zero-overpressure  radii  as  functions  of  time  based  on 
fits  to  analytical  solutions.  These  fimctions  may  be  inverted  using  Newton’s  method  to  obtain  the  arrival 
times.  We  developed  alternative  fits  to  experimental  data  for  a  somewhat  different  standard  which 
provides  the  shock  arrival  time  and  positive  phase  duration  as  functions  of  position.  These  can  be  used 
directly  to  compute  shock  and  zero-overpressure  arrival  times. 


2 


The  ground  plane  is  assumed  to  be  perfectly  flat  and  rigid  and  to  effectively  double  the  explosive 
weight  of  a  stack.  The  latter  is  also  a  worst-case  assumption.  Waves  from  multiple  sources  are  assumed 
to  propagate  independently  from  one  another,  and  coalescence  is  assumed  to  occur  wherever  and  whenever 
the  shock  associated  with  one  blast  wave  encroaches  into  the  positive  overpressure  phase  associated  with 
another  blast  wave. 


Under  these  assumptions,  regardless  of  the  spatial  and  temporal  separation  between  the  sources,  at 
sufficiently  large  distances  from  the  charges,  coalescence  will  be  predicted.  Blast  waves  from  actual 
detonating  stack  arrangements  do  not  exhibit  this  property  because  they  do  not  propagate  independently. 
Rather,  portions  of  some  of  them  propagate  in  air  processed  by  their  predecessors.  An  overtaking  wave 
may  first  encounter  the  negative  overpressure  phase  of  the  wave  being  overtaken  causing  it  to  decelerate. 
The  overtaking  wave  must  be  somewhat  stronger  than  the  wave  being  overtaken  or  it  will  not  be  able  to 
penetrate  the  negative  overpressure  phase  and  coalescence  cannot  occur.  Once  this  phase  has  been 
penetrated  and  the  overtaking  wave  has  entered  the  positive  phase  of  the  wave  being  overtaken, 
coalescence  is  assured  but  occurs  only  after  some  additional  propagation.  For  these  reasons,  the  BWACO 
coalescence  algorithm  can  be  expected  to  predict  coalescence  where,  in  actuality,  it  does  not  occur.  A 
margin  of  safety  is  provided  by  prediction  of  a  larger  hazard  area  than  actually  exists. 


We  considered  several  methods  of  combining  peak  overpressures  from  N  coalesced  waves.  These 
include  simple  superposition, 

N 

Ap=ZApi, 

i=l 


the  full  LAMB  algorithm, 
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i=l 
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or  use  of  the  pressure  produced  by  the  total  explosive  weight  of  all  stacks  associated  with  the  coalesced 
waves  combined  at  their  center  of  charge.  This  last  approach  is  described  in  greater  detail  in  section  5.5. 

4.  USING  BWACO 

In  order  to  run  BWACO,  the  user  supplies  a  description  of  the  distributed  ammxinition  stacks  and,  if 
desired,  of  the  region  of  interest.  The  program  divides  the  region  into  discrete  stations.  It  is  recognized 
that,  generally,  any  stack  may  act  as  the  initial  d(Mior  and  it  is  necessary  to  determine  the  worst-case 
loading  (highest  combined  peak  overpressure)  at  each  station  considering  aU  initial  donors.  For  each 
initial  donor,  BWACO  determines  the  timing  of  detonation  of  each  of  the  other  stacks.  It  then  applies  the 
coalescence  criterion  and  computes  both  the  total  equivalent  weight  of  the  explosive  in  the  stacks 
contributing  to  the  blast  and  the  peak  overpressure  at  each  station.  It  retains  the  values  associated  with 
the  highest  peak  overpressure  for  any  initial  donor.  The  results  are  reported  in  graphical  form. 

User  specifications  are  supplied  by  means  of  an  appropriately  formatted  input  file.  Each  ammunition 
stack  is  defined  by  specifying  its  coordinates  in  a  user-defined  Cartesian  system,  its  net  equivalent 
explosive  weight,  a  propagation  velocity  for  communication  of  mass  detonation  to  neighboring  stacks,  and 
whether  or  not  the  stack  may  act  as  an  initial  donor.  Any  convenient  Cartesian  coordinate  system  may 
be  used  with  distances  specified  in  meters.  The  explosive  weight  for  the  TNT  equivalent,  in  pounds,  of 
all  the  explosive  in  the  stack  is  required,  as  is  the  propagation  velocity  in  m/s.  Usually,  in  order  to  allow 
detennination  of  the  worst-case  blast  loading,  all  of  the  stacks  should  be  allowed  to  act  as  initial  donors. 
However,  this  input  option  provides  the  flexibility  required  to  address  special  problems  such  as  simulating 
experiments  with  known  donors.  The  extent  of  the  region  of  interest  is  specified  using  the  same 
coordinate  system  defined  for  the  stacks.  If  no  region  is  specified,  BWACO  selects  a  region  limited  to 
the  area  in  which  the  total  equivalent  weight  of  the  explosive  in  all  the  stacks  produces  peak  blast 
overpressures  above  the  allowable  value  for  inhabited  buildings  (unless  similar  considerations  for  the 
individual  stacks  dictate  a  larger  region).  The  program  then  positions  3,600  stations  within  the  region. 

BWACO  cycles  through  all  possible  initial  donors  as  specified  by  the  user.  For  each  initial  donor, 
using  the  given  propagation  velocities,  it  computes  the  times  at  which  the  other  stacks  detonate.  It  applies 
the  coalescence  criterion  and,  if  necessary,  combines  the  peak  overpressures  at  each  station. 

BWACO  provides  different  graphical  representations  of  the  results  depending  on  the  version  used. 
The  coalescence  map  is  a  plot  showing  areas  of  coalescence  within  the  region  of  interest  and  indicating. 
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in  the  Cray  version,  the  number  of  waves  which  have  coalesced  or,  in  the  PC  version,  the  total  equivalent 
explosive  weight  contributing  to  the  blast.  The  peak  overpressure  map,  available  only  in  the  PC  version, 
is  a  contour  plot  showing  the  worst-case  peak  overpressure  from  the  blast  in  the  region  of  interest.  The 
contour  intervals  correspond  to  the  five  permissible  exposure  levels  defined  in  AR  385-64. 

5.  EVOLUTION  OF  THE  BWACO  MODEL 

5.  1  Review  of  Experimental  Results.  Zaker  (1969)  reported  results  of  an  extensive  analytical  and 
experimental  study  of  the  coalescence  of  blast  waves  produced  by  pairs  of  sequentially  detonating 
Composition  C-4  charges  having  a  total  weight  of  two  pounds.  In  the  analysis,  the  charges  were  assumed 
to  be  located  at  the  same  point,  while  in  the  experiments  they  were  separated  by  10  inches  and  a  steel 
barrier  to  prevent  sympathetic  detonation.  Pressure  was  measured  along  a  lateral  line  equidistant  from  the 
charges  and  along  an  axial  line.  The  analysis  was  assumed  to  be  applicable  to  the  experiments  along  the 
lateral  line.  The  axial  values  measured  depend  on  the  separation  between  the  charges  and  are  affected  by 
the  presence  of  the  barrier,  while  the  lateral  values  are  approximately  independent  of  the  separation  and 
are  predicted  by  the  analysis. 

Zaker  considered  charge  weight  ratios  for  successively  detonated  charges  of  1:2,  1:1,  and  2:1  and 
nominal  initiation  delay  times  ranging  from  0.8  to  5.7  ms  corresponding  to  0.60  to  4.11  ms/lb^^  using 
conventional  blast  scaling.  From  this  it  can  be  inferred  that  an  equivalent  TNT  weight  of  about  2.62  lb 
was  used  in  the  scaling.  The  blast  environment  was  monitored  using  pressure  gauges  on  the  axis  of 
charge  centers  and  in  the  lateral  direction  peipendicular  to  this  axis  at  the  center  of  charge  to  a  distance 
of  aRjroximately  58  ft  (42  ft/lb^^).  This  distance  is  a  little  shorter  than  that  at  which  the  overpressure 
from  the  combined  charges  decays  to  about  0.9  psi.  This  represents  the  smallest  permissible  exposure 
level  (which  applies  to  inhabited  buildings)  defined  in  AR  385-64. 

Pressure  records  showing  the  coalescence  process  were  presented  by  Zaker.  These  indicate  that 
substantial  propagation,  covering  most  of  the  field  of  observation,  may  occur  as  coalescence  progresses. 
He  observed  that  the  peak  overpressures  associated  with  coalesced  waves  were  essentially  the  same  as 
those  produced  by  single  charges  of  the  same  total  explosive  weight. 

For  two  charges  of  equal  mass,  Zaker  found  that  a  "tendency"  to  coalescence  in  the  lateral  direction 
occurs  for  delays  of  less  than  4,3  ms  (3.2  ms/lb’^).  The  term  "tendency"  was  not  clearly  defined  but  it 
may  be  assumed  to  mean  observation  of  decreasing  intervals  between  wave  peaks  with  increasing  distance 
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from  the  charge  center.  This  is  actually  very  similar  to  the  BWACO  criterion.  He  observed  a  tendency 
to  coalescence  in  the  axial  direction  with  all  of  the  delays  considered  (up  to  5.7  ms  or  4.1  ms/lb^^).  By 
comparing  lateral  and  axial  data,  he  estimated  that  the  effect  of  the  charge  separation  and  the  barrier  on 
the  interval  between  peaks  in  the  axial  direction  is  equivalent  to  an  additional  delay  of  1.8  ms 
(1.3  ms/lb^^).  Thus,  axial  coalescence  may  be  assumed  to  occur  for  delays  less  than  6.1  ms 
(4.5  msAb^^).  This  result  may  include  significant  effects  of  the  barrier. 

For  charges  of  unequal  mass,  Zaker  observed  that  coalescence  persists  at  longer  delays  with  a  weight 
ratio  of  1:2  and  vanishes  at  shorter  delays  with  a  weight  ratio  of  2:1  compared  to  equal  charges.  (Thus, 
the  criterion  given  in  AR  385-64  uses  more  conservative  values  of  4.0  and  5.6  msAb^^,  respectively.) 
7^tfpr  also  performed  some  experiments  with  three  equal  charges.  He  observed  that  the  third  pulse  tends 
to  overtake  the  second  before  the  second  overtakes  the  first. 

5.2  Simulation  of  the  Experiments  Using  the  1-KT  Standard.  In  order  to  simulate  Zaker’s 
equal-charge  experiments,  we  made  BWACO  computations  for  two  1.31 -lb  "stacks"  separated  by  0.254  m, 
specifying  only  one  possible  initial  donor.  Specific  delay  times  were  obtained  using  appropriate  values 
of  the  propagation  velocity.  The  region  of  interest  was  specified  as  extending  to  18.0  m  from  the  center 
of  charge  in  each  direction,  corresponding  to  the  region  of  experimental  observation.  The  delay  times 
associated  with  the  vanishing  of  coalescence  at  the  edges  of  the  region  may  be  compared  with 
experimental  values. 

Coalescence  maps  of  regions  of  the  surface  for  various  delay  times  are  plotted  in  Figure  1.  The  donor 
(solid  circle)  and  acceptor  (solid  square)  are  plotted  near  the  center  of  the  map  and  overlap.  Regions  of 
coalescence  are  shaded  with  dots.  Unshaded  areas  experience  two  independent  waves,  while  shaded  areas 
experience  a  single  coalesced  wave.  For  delays  less  than  3.3  ms,  coalescence  within  the  region  of  interest 
is  predicted  in  all  directions.  For  delays  greater  than  4.0  ms  coalescence  in  the  lateral  (x)  direction 
vanishes.  Thus,  with  a  delay  of  4.3  ms  (Zaker’s  limit  for  lateral  coalescence),  coalescence  in  the  lateral 
direction  is  not  predicted.  Coalescence  in  the  axial  (y)  direction  persists  until  the  delay  exceeds  4.7  ms. 
This  may  be  compared  to  Zaker’s  value  of  6.1  ms.  The  difference  between  the  delays  for  which 
coalescence  vanishes  in  the  axial  and  lateral  directions  respectively  is  only  0.7  ms,  which  is  much  less  than 
Zaker’s  value  of  1.8  ms. 

In  each  direction,  vanishing  of  coalescence  is  predicted  at  shorter  delays  than  observed  in  or  estimated 
from  the  Zaker  experiments.  The  discrepancy  is  most  significant  in  the  axial  direction.  This  result,  which 
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is  unacceptable,  contradicts  our  expectations  for  the  bias  in  BWACO’s  coalescence  algorithm  and 
suggested  to  us  that  the  1-KT  Standard  might  be  inaccurate. 

5.3  Alternatives  to  the  1-KT  Standard  Fits.  Kingery  and  coworicers  have  reported  air  blast 
parameters,  including  shock  arrival  time,  positive  phase  duration  and  peak  overpressure,  versus  distance 
for  large-scale  hemispherical  TNT  surface  bursts  (Kingery,  Keefer,  and  Day  1962;  Kingery  and  Pannill 
1964;  and  Kingery  1966).  This  data  provides  a  basis  for  evaluating  the  applicability  of  the  1-KT  Standard. 
Comparison  with  the  1-KT  Standard  shock  arrival  time  is  shown  in  Figure  2.  There  is  close  agreement 
except  in  a  tiny  region  very  near  the  source.  Comparison  with  the  1-KT  Standard  positive  phase  duration 
is  shown  in  Figure  3.  The  values  from  the  standard  are  generally  too  short  at  distances  of  interest  This 
would  tend  to  suppress  the  prediction  of  coalescence,  producing  the  results  we  observed.  Comparison  with 
the  1-KT  Standard  peak  overpressure  is  shown  in  Figure  4.  Zaker’s  small-scale  results  are  also  plotted. 
The  1-KT  Standard  pressures  are  generally  lower  than  the  measured  values. 

Alternative  fits  which  give  shock  arrival  time  and  positive  phase  duration  as  functions  of  distance  were 
developed.  These  retain  the  appropriate  asymptotic  behavior  at  large  distances  from  the  source.  The  fit 
for  shock  arrival  time  is  of  the  form 


tj  =  ag[r-batanh(_L)], 

» 


where  a.,  b„  and  r„  are  constants.  The  fit  for  positive  phase  duration  is  of  the  form 

a*  a  A 

At^=a,(l--^)[ln(r2+r5]^^. 


where  a^.,  b^.,  c^.  and  r.^  are  constants.  This  fit  is  not  used  for  very  small  values  of  r,  for  which  the  value 
of  the  positive  phase  duration  is  held  constant.  Kingery  and  Pannill  (1964)  provided  a  fit  for  peak 
overpressure.  These  fits  are  also  plotted  in  Figures  2  through  4  and  provide  improved  agreement  with 
measurements  in  each  case.^ 


t 


After  implementation 


of  this  standard,  we  learned  of  the  existence  of  a  TNT  standard. 
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Figure  2.  Comparison  of  the  1-KT  standard  with  the  alternative  standard  for  shock  arrival  time. 
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Figures.  Comparison  of  the  1-KT  standard  with  the  alternative  standard  for  positive  phase  duration. 
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5.4  Simulation  of  the  Experiments  Using  the  Alternative  Standard.  Coalescence  maps  obtained  using 
the  alternative  standard  are  shown  in  Figure  5.  Lateral  coalescence  vanishes  for  delays  exceeding  about 

5.1  ms,  which  is  significantly  greater  than  Taker’s  observation  of  4.3  ms.  In  the  axial  direction, 
coalescence  vanishes  for  delays  of  5.9  ms  or  greater.  This  is  only  slightly  less  than  Taker’s  estimate  of 

6.1  ms.  There  is  no  significant  change  in  the  difference  between  the  axial  and  radial  values.  It  is  still 
less  than  half  that  determined  by  Taker.  The  predictions  obtained  with  the  alternative  standard  do  not  miss 
significant  regions  of  coalescence  and  are  acceptable. 

As  we  have  noted,  at  sufficiently  large  distances  from  the  stacks,  coalescence  will  be  predicted.  For 
the  present  conditions,  this  becomes  apparent  upon  expanding  the  regions  of  interest  while  maintaining 
the  other  parameters  as  shown  in  Figure  6.  Although  coalescence  is  identified  at  large  distances,  the 
pressures  produced  are  too  low  to  be  of  concern.  BWACO  was  configured  to  suppress  identification  of 
coalescence  where  the  combined  peak  overpressure  is  below  0.9  psi  (the  maximum  permissible  exposure 
level  for  inhabited  buildings  subjected  to  blast  waves  from  large  charges).  The  effect  on  the  coalescence 
maps  is  to  eliminate  shading  where  distance  from  the  charges  is  sufficiently  large  as  illustrated  in  Figure 
7.  It  appears  that  the  pressures  combined  using  simple  superposition  are  somewhat  high,  as  the  extent  of 
the  region  of  significant  pressure  is  greater  than  that  associated  with  a  single  charge  of  the  same  total 
weight.  As  a  result  of  the  expanded  extent,  coalescence  is  observed  at  somewhat  longer  delays  than  those 
associated  with  the  experimental  region  of  significant  pressure  specified  in  the  earlier  computations.  (See 
Figure  5.) 

5.5  Effect  of  Pressure  Combination  Algorithms.  To  assess  the  effect  of  the  pressure  combination 
algorithms  on  the  extent  of  the  region  of  significant  pressure,  we  made  computations  with  from  2  to  10 
stacks,  located  at  the  same  point,  having  the  same  total  equivalent  explosive  weight  (2.62  lb).  The  ratio 
of  the  radius  of  the  region  of  significant  pressure  produced  by  the  detonation  of  n  stacks  to  that  produced 
by  the  detonation  of  one  stack  is  plotted  as  a  function  of  n  in  Figure  8.  Contrary  to  Taker’s  observation, 
using  simple  superposition,  we  obtain  substantial  increase  in  the  extent  of  the  region  of  significant  pressure 
with  increasing  number  of  stacks  even  though  the  total  explosive  weight  remains  constant. 

The  results  obtained  with  the  LAMB  algorithm  (also  shown  in  Figure  8)  are  virtually  identical  to  those 
obtained  using  simple  superposition.  The  LAMB  algorithm  has  been  shown  to  produce  correct  results  for 
the  head-on  collision  of  two  equal  blast  waves  and  overestimates  the  pressure  in  other  cases  (Hikida  and 
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Figure  8.  Comparison  of  the  effects  of  the  pressure  combination  algorithm  on  the  region  of  significant 
pressure  as  a  function  of  the  number  of  charges. 

Needham  1981;  Erode  1977).  In  the  case  of  overtaking  waves,  the  magnitude  of  the  overestimate  appears 
significant. 

As  an  alternative  to  the  simple  supeiposition  and  LAMB  algorithms,  we  computed  the  combined 
pressure  by  first  combining  the  stacks  which  produce  the  coalesced  wave  into  a  single  stack  positioned 
at  the  center  of  charge.  This  aj^roach,  suggested  by  Taker’s  observation,  is  applicable  where  the  angles 
between  coalescing  waves  are  small,  as  occurs  at  any  point  whose  distance  from  the  stacks  significantly 
exceeds  the  distance  between  stacks.  It  produces  regions  of  significant  overpressure  which  do  not  increase 
with  the  number  of  stacks.  This  result  is  also  shown  in  Figure  8. 
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5.6  Simulation  of  the  Experiments  Using  the  Alternative  Standard  and  Pressure  Combination 


Algorithm.  The  way  in  which  the  alternative  pressure  combination  algorithm  modifies  the  results  obtained 
with  two  equal  charges  is  shown  in  Figure  9.  As  expected,  the  region  of  significant  pressure  is  reduced 
in  comparison  to  that  shown  in  Figure  7.  The  delays  at  which  coalescence  vanishes  are  slightly  increased 
compared  with  those  associated  with  the  analysis  shown  in  Figure  5  because  the  region  of  significant 
pressure  is  larger  than  the  region  of  experimental  observation. 

We  also  made  computations  with  unequal  charges  in  ratios  of  1:2  and  2:1  with  a  constant  total 
equivalent  weight  of  2.62  lb.  Figure  10  shows  coalescence  maps  for  the  1:2  case  at  delays  ranging  from 
4.0  to  6.5  ms.  There  is  still  a  small  region  of  coalescence  at  the  longest  delay.  Figure  11  shows 
coalescence  maps  for  the  2:1  case  at  delays  ranging  from  3.0  to  5.5  ms.  In  this  case,  coalescence  has 
completely  vanished  at  the  longest  delay.  Figure  12  compares  results  from  each  of  these  cases  with  the 
equal-charge  case  at  two  different  delays.  This  shows  that  the  persistence  of  coalescence  decreases  as  the 
ratio  of  the  donor  weight  to  the  acceptor  weight  increases.  The  result  is  consistent  with  Zaker’s 
observations. 

Figure  13  shows  coalescence  maps  obtained  in  a  simulation  of  Zaker’s  three-charge  experiment  The 
size  of  each  shading  dot  is  proportional  to  the  number  of  waves  coalesced  at  the  corresponding  station. 
However,  it  cannot  be  ascertained  which  waves  have  coalesced.  For  short  delays,  coalescence  of  all  three 
waves  is  predicted  in  all  directions,  although  two-wave  coalescence  occurs  first  (closer  to  the  center  of 
charge).  As  the  delay  is  increased,  three-wave  coalescence  begins  to  vanish,  starting  in  the  negative  axial 
direction.  Two-wave  coalescence  persists  for  longer  delays.  The  results  are  consistent  with  Zaker’s 
observation  that  the  third  pulse  tends  to  overtake  the  second  before  the  second  overtakes  the  first 

5.7  Results  With  Multiple  Initial  Donors.  To  this  point,  all  the  computations  have  specified  a  single 
donor.  Each  resultant  coalescence  map  shows  the  coalescence  pattern  associated  with  a  single  chain  of 
events.  However,  BWACO  is  generally  intended  to  be  used  with  all  stacks  specified  as  possible  initial 
donors.  The  coalescence  map  produced  in  this  case  is  a  composite  which  reflects  the  worst  case  loading 
at  each  point  in  the  region  of  interest  The  effect  of  allowing  each  stack  to  act  as  an  initial  donor  in  the 
Zaker  simulations  was  examined  in  order  to  check  the  code’s  multiple  donor  algorithm. 
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The  coalescence  patterns  produced  in  simulations  of  two  equal  charges  are  shown  in  Figure  14.  As 
expected,  symmetry  is  produced,  with  the  single-donor  coalescence  regions  (shown  in  Figure  9)  reflected 
across  the  horizontal  axis. 

The  coalescence  regions  produced  with  two  unequal  charges  are  not  symmetrical  about  the  horizontal 
axis  as  shown  in  Figure  15.  These  are  combinations  of  the  single-donor  regions  (shown  in  Figure  10)  with 
reflections  of  the  single-donor  regions  (shown  in  Figure  11)  about  the  horizontal  axis. 

Symmetry  about  the  horizontal  axis  also  results  with  three  equal  charges,  as  shown  in  Figure  16.  In 
this  case,  however,  the  region  of  coalescence  is  greater  than  that  associated  with  reflection  of  the  single 
donor  regions  (shown  in  Figure  13)  because  coalescence  produced  by  the  center-charge  donor  dominates. 
In  fact,  coalescence  of  the  two  waves  from  the  outer  charges  persists  at  all  delays  since  they  detonate 
simultaneously  when  the  central  charge  is  the  donor. 

6.  PERSONAL  COMPUTER  VERSION  GRAPHICS 

When  BWACO  was  adapted  for  use  on  the  personal  computer,  the  coalescence  map  was  enhanced 
with  the  addition  of  contours  of  the  total  explosive  weight  associated  with  the  coalesced  waves  and  a 
capability  for  producing  peak  overpressure  contour  plots  was  added.  Examples  of  these  maps  are  shown 
in  Figure  17.  The  maps  may  be  obtained  either  on  the  PC  display  or  in  hard-copy  form. 

Regions  in  which  coalescence  is  detected  are  hatched  in  the  coalescence  map.  Each  contour  is  labeled 
to  indicate  the  total  explosive  weight  contributing  to  the  worst-case  blast  environment  within  that  contour. 
However,  the  outermost  contour  line  is  further  limited  by  the  overpressure  limit  for  inhabited  buildings. 
An  inhabited  building  distance  circle  centered  at  the  center  of  the  charge  is  also  shown.  The  legend  lists 
all  the  contour  levels  present  in  the  plot.  The  highest  contour  level  represents  the  maximum  effective 
explosive  weight  for  the  configuration. 

In  the  peak  overpressure  map,  contours  are  labeled  to  indicate  the  worst-case  peak  oveipressure  at  each 
point  within  the  region  of  interest.  The  contour  levels  are  those  specified  in  AR  385-64.  The  lowest  level 
varies  between  0.9  and  1.2  psi,  depending  on  the  total  explosive  weight.  The  limits  of  the  blast  hazard 
area  determined  in  accordance  with  the  regulations  are  also  shown.  The  legend  lists  the  stack  positions 


24 


25 


Distance  (m)  Distance  (m)  _ _ _ 


28 


and  weights,  along  with  center  of  charge  and  total  weight,  compares  the  regulatory  blast  hazard  with  that 
determined  by  BWACO,  and  lists  the  contour  levels. 


7.  SIMULATIONS  OF  REPRESENTATIVE  LARGE-SCALE  CONHGURATIONS 

In  the  large-scale  scenarios  that  BWACO  is  intended  to  treat,  initiation  delay  depends  on  the 
separation  between  stacks.  The  simulations  described  in  the  following  paragraphs  were  intended  to 
illustrate  the  effects  of  the  separation  between  stacks  and  the  distribution  of  explosive  weight  among  stacks 
when  all  stacks  act  as  initial  donors.  Most  of  them  were  made  using  a  total  equivalent  weight  of 
100,000  lb  distributed  over  two  or  more  stacks.  For  this  total  weight,  regulations  permit  an  inhabited 
building  exposure  of  1.2  psi  or  less.  The  velocity  of  propagation  between  stacks  was  fixed  at  2,500  m/s. 

Results  obtained  with  two  stacks  of  50,000  lb  each  are  shown  in  Figures  18  and  19.  At  the  smallest 
separation  distance  (50  m),  the  regulatory  blast  hazard  region  nearly  coincides  with  the  BWACO 
prediction.  As  the  separation  is  increased,  significant  regions  in  which  the  pressure  is  below  1.2  psi 
appear.  In  these  cases,  requiring  consideration  of  the  blast  produced  by  the  total  explosive  weight  of  both 
stacks  is  clearly  too  restrictive.  At  the  larger  separation  distances,  regions  of  coalescence  persist  only  near 
the  horizontal  symmetry  axis.  Here,  consideration  of  the  total  weight  is  necessary.  The  variation  of  the 
blast  hazard  areas  with  separation  is  summarized  in  Table  1.  The  regulatory  hazard  area  increases  with 
separation  while  that  predicted  by  BWACO  is  minimized  somewhere  between  100  and  150  m. 


Table  1.  Variation  of  Blast  Hazard  Areas  with  Separation  for  Two  50,000-lb  Stacks 


Blast  Hazard  Area 

Separation 

AR  385-64 

BWACO 

(m) 

(km^) 

(km^) 

50 

1.06 

0.97 

100 

1.12 

0.86 

150 

1.18 

0.86 

200 

1.23 

0.88 
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Figure  18.  Coalescence  maps  for  large-scale  simulations  with  two  50,000-lb  stacks. 
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We  also  simulated  the  effects  of  varying  the  distribution  of  explosive  weight  betweai  two  stacks 
sq)arated  by  150  m.  The  results,  shown  in  Figures  20  and  21,  indicate  very  little  diange  in  the 
coalescaice  region  along  the  horizontal  axis  as  tte  weight  is  redistributed.  The  hazard  area  increases 
subtly  for  increasingly  unequal  weight  distributions  as  shown  in  Table  2. 


Table  2.  Variation  of  Blast  Hazard  Areas  widi  Explosive  Weight  Distribution  for  Two  Stadts 
Having  a  Total  Weight  of  100,000  lb  Separated  by  150  m 


Weight  Ratio 
(lb:lb) 

Blast  Hazard  Area 

AR  385-64 
(km^) 

BWACO 

(km^) 

50,000  :  50,000 

1.18 

0.86 

40,000  ;  60,000 

1.18 

0.85 

30,000  :  70,000 

1.18 

0.85 

20,000  :  80,000 

1.18 

0.89 

10,000  :  90,000 

1.18 

0.93 

The  results  observed  with  four  25,000-lb  stacks  are  similar  to  those  seen  with  two  50,000-lb  stacks. 
They  are  shown  in  Figures  22  and  23.  Coalescence  of  waves  from  all  four  stacks  and  from  three  of  the 
four  is  only  exhibited  at  the  50-m  separation  distance.  These  regions  ^pear  as  small  islands  lying 
beyond  die  general  regions  of  coalescence.  At  the  100-m  separation,  only  two-stack  coalescence  results. 
As  the  separation  is  further  increased,  this  coalescence  becomes  limited  to  the  vicinity  of  the  vertical  and 
horizontal  symmetry  axes.  The  effect  of  separation  on  the  blast  hazard  areas  is  given  in  Table  3.  Again, 
BWACO  predicts  a  minimum  hazard  near  the  100-m  separation. 


Table  3.  Variatiwi  of  Blast  Hazard  Areas  with  Separation  for  Four  25,000-lb  Stacks 


Blast  Hazard  Area 

Separation 

AR  385-64 

BWACO 

(m) 

(km^) 

(km^) 

50 

1.13 

0.85 

100 

1.25 

0.71 

150 

1.37 

0.73 

The  effect  of  the  number  of  stacks  was  also  investigated.  If  each  stack  in  the  4-stack  arrangemait  is 
divided  into  4  equal  stacks  (6,250  lb  ea)  which  are  positioned  such  that  dieir  center  of  charge  lies  at  the 
position  of  the  original  stack  while  retaining  a  uniform  separation,  an  arrangement  of  16  stacks  with  half 
the  original  separation  results.  The  coalescence  and  peak  overpressure  maps  produced  are  shown  in 
Figures  24  and  25.  Each  map  may  be  compared  lo  the  corresponding  map  in  Figures  22  or  23.  A 
significant  advantage  to  smaller  subdivisions  is  indicated.  The  effect  of  separation  on  the  maximum 
effective  explosive  weight  and  the  blast  hazard  areas  is  given  in  Table  4.  -  Comparison  with  Table  3  shows 
a  reduction  in  the  predicted  blast  hazard  area  but  an  increase  in  die  regulatory  area.  The  minimum 
predicted  hazard  area  occurs  near  the  SO-m  separation. 


Table  4.  Variation  of  Maximum  Effective  Explosive  Weight  and  Blast  Hazard  Areas 
ivith  Separation  for  Sixteen  6,250-lb  Stacks 


Separation 

(m) 

Maximum  Effective 
Explosive  Weight 

Ob) 

Blast  Hazard  Area 

AR  385-64 
(km^) 

BWACO 

(km^) 

25 

56,250 

1.18 

0.64 

50 

25,000 

1.37 

0.50 

75 

25,000 

1.57 

0.56 

Simulations  of  representative  storage  configurations  with  6  and  10  unequal  stacks  demonstrate  the 
utility  of  BWACO.  In  each  case,  the  stacks  have  TNT  equivalent  weights  between  20,000  and  30,(XX)  lb 
and  are  separated  from  their  neighbors  by  a  distance  of  100  m. 

The  six  stacks  have  a  total  weight  of  150,000  lb.  The  coalescence  and  peak  overpressure  m^  are 
shown  in  Figure  26.  The  weight  contours  in  the  coalescence  map  are  hard  to  distinguish  but  their  levels 
are  listed  in  the  legend.  The  highest  weight  effective  within  the  regicni  of  significant  {xessure  is 
99,000  lb.  This  is  significantly  less  than  the  total  weight.  The  level  of  tte  outermost  contour  shown  on 
the  peak  overpressure  map  is  1.1  psi.  This  corresponds  to  the  regulatory  lequiremoit  for  the  150,000-lb 
total  weight  The  blast  hazard  area  predicted  by  BWACO  using  tiiis  pressure  level  is  less  than  half  drat 
specified  by  the  regulation.  If  the  larger  peak  overpressure  level  corresponding  to  the  maximum  effective 
weight  (99,000  lb)  were  used,  an  even  more  favorable  comparison  would  result 
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The  ten  stacks  have  a  total  weight  of  250,000  lb,  requiring  consideration  of  peak  overpressures  down 
to  0.9  psi.  The  coalescence  and  peak  overpressure  maps  are  shown  in  Figure  27. 

The  highest  weight  effective  within  the  region  of  significant  pressure  is  134,000  lb.  The  ratio  of  the 
BWACO  hazard  area  to  the  regulatory  hazard  area  is  even  smaller  than  in  the  six-stack  case.  Again,  this 
comparison  would  become  still  more  favorable  if  the  larger  overpressure  corresponding  to  the  maximum 
effective  weight  were  used. 

8.  SUMMARY  AND  CONCLUSIONS 

In  this  report,  we  have  explained  the  assumptions  underlying  the  BWACO  algorithms,  documented 
the  evolution  of  BWACO  based  on  comparisons  with  available  experimental  data,  and  demonstrated  the 
application  of  BWACO  to  typical  large-scale  ammunition  storage  configurations. 

Assumptions  were  made  in  order  to  determine  the  order  and  timing  of  the  detonation  of  the  stacks 
following  the  detonation  of  any  initial  donor  stack,  to  establish  a  criterion  for  detecting  coalescence,  and 
to  determine  the  combined  pressure  associated  with  a  number  of  coalesced  waves  at  a  point. 

Comparison  of  preliminary  results  with  experimental  data  obtained  by  Zaker  led  to  replacement  of  the 
standards  initially  used  for  the  description  of  blast  waves.  Use  of  the  1-KT  Standard  was  found  to 
produce  predictions  which  tended  to  miss  detection  of  coalescence  where  the  experiments  showed  that 
coalescence  occurred.  It  was  replaced  with  a  standard  based  on  experimental  data  reported  by  Kingery. 
In  addition,  a  pressure  combination  algorithm  based  on  Zaker’s  observation  was  foimd  to  produce  more 
satisfactory  results  than  other  algorithms. 

BWACO  has  been  adapted  for  the  personal  computer  with  enhanced  graphical  representations.  As 
currently  configured,  BWACO  provides  a  means  of  assessing  the  blast  environment  associated  with  the 
sequential  detonation  of  an  arbitrary  arrangement  of  ammunition  stacks.  The  limitations  imposed  by  the 
assumptions  have  not  been  assessed  in  realistic  configurations. 

Applications  to  a  number  of  problems  representative  of  typical  ammunition  storage  configurations  were 
detailed.  The  results  indicate  that  regions  of  significant  pressure  associated  with  the  coalescence  of  blast 
waves  fiom  distributed  ammunition  stacks  may  be  less  extensive  than  corresponding  regions  associated 
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Figure  27.  Coalescence  and  peak  overpressure  maos  for  a  large-scale  simulation  with  ten  unequal 
stacks  having  a  total  weight  of  250,000  lb. 
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with  regulatory  requirements.  An  advantage  associated  with  the  distribution  of  ammunition  into  smaller 
subdivisions  was  also  demonstrated. 

If  it  proves  advantageous,  opportunities  to  apply  more  sophisticated  models  can  be  exploited.  Areas 
where  improvement  is  possible  include  the  propagation  model  and  the  pressure  combination  algorithm, 
as  well  as  the  addition  of  terrain  effects  and  fragment  hazard  predictions. 
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APPENDIX  A: 

DETAILED  DESCRIPTION  OF  BWACO  SOURCE  CODE 
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BACKGROUND 


The  BWACO  program  was  writtoi  to  {xedict  regions  of  blast  wave  ooalescoice  associated  with  the 
sequential  mass  detonation  of  distributed  ammunition  stadns  positkHied  as  specified  by  the  user.  It 
accomplishes  this  by  assuming  that  blast  waves  {vopagate  indqjendoifiy  of  one  anotl^r  and  that 
coalescence  occurs  wherever  and  whenever  the  shock  associated  with  one  blast  wave  encroaches  into  the 
positive  overpiessuie  i^iase  of  another  blast  wave.  It  further  assumes  that  communication  of  mass 
detonation  between  stacks  is  omtrolled  by  a  user-supplied  propagation  vdodty  (whidi  mi^  represent 
the  velocity  of  fiagmraits  from  the  donor  stack).  It  uses  this  {Hopagatkm  vdodty  to  determine  the  order 
and  timing  of  detonation  of  the  stacks.  It  recognizes  that  die  identity  of  die  initial  donor  is  not  g^rally 
known  and,  unless  a  single  donor  has  been  specified  by  the  user,  considers  all  possible  donors  in  order 
to  predict  the  worst-case  environment  at  any  point.  It  reports  the  peak  overpressure  and  effective 
explosive  weight  ctmtrihuting  to  this  environment  throughout  the  region  of  interest. 

Although  the  vast  majority  of  BWACO  represents  new  code,  a  few  of  its  FORTRAN  subroutines  and 
functions  were  adapted  from  the  DNA  Nuclear  Blast  Standard  (1  KT)  (hereinafter  referred  to  as  die  1-KT 
Standard),  which  describes  the  blast  environment  produced  by  a  1-KT  nuclear  explosion,  and  from  the 
Low  Altitude  Multiple  Burst  (LAMB)  model,  which  describes  the  Wast  raivirwiment  produced  by  multiple 
nuclear  explosions  of  artHtrary  yield.  As  a  result,  mudh  of  the  terminology  used  to  discuss  the  program 
comes  from  the  nuclear  blast  vocabulary.  For  example,  the  net  equivalent  explosive  weight  in  an 
ammunition  stack  is  referred  to  as  its  yield  and  detonation  of  a  stadc  is  called  a  burst. 

The  organization  of  BWACO’s  subroutines  and  functions  is  diagrammed  in  Hgure  A-1  (where 
remaining  adaptations  from  the  1-KT  Standard  or  LAMB  are  rmdered  white  on  black).  The  routines  ate 
divided  into  six  modules.  The  interface  module  directs  BWACO’s  major  activities,  configuring  the  user 
session,  setting  the  constants,  and  interfacing  with  the  user  by  providing  a  mmu  of  options.  These  include 
problran  setup  and  execution,  plotting,  and  recorrfiguration  of  the  session.  The  executicm  module  reads 
the  user’s  irqrut,  determines  the  region  of  interest,  distributes  stations  within  it,  and  executes  the 
cnmpitatinn.  The  plotting  module  interfaces  with  the  user  to  produce  plots  tm  request  The  propagaticxi 
module  determines  the  order  and  timing  of  detonatkm  propagation  ditough  tire  stadts.  The  coalescence 
module  determines  the  waveform  variables  and  ^^lies  the  coalescence  criterion  at  each  staticm.  The 
waveform  module  computes  the  requited  scale  factors  and  determines  the  requited  waveform  parameters 
fiom  the  jqrplicable  standard. 
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The  following  discussion  details  most  of  the  routines,  describing  the  program  from  the  inside  out. 
B WACO’S  logic  is  difficult  to  understand  and  explain  and  the  inside-out  i^proach  facilitates  exposition 
of  its  operation. 

WAVEFORM  MODULE 

These  routines  use  the  applicable  blast  standard  to  describe  the  diock  wave,  die  boundary  between  the 
positive  and  negative  phases  of  die  blast  wave  (point  of  zero  overpressure)  and  the  peak  values  of 
overpressure,  density,  and  particle  velocity  associated  with  a  burst. 

Subroutine  at62,  adapted  from  the  1-KT  Standard,  is  shown  in  Figure  A-2.  It  computes  ambient 
pressure  (Pa),  sound  speed  (cm/s),  density  (g/cm^)  and  temperature  (K)  from  the  1962  U.S.  standard 
temperature  atmosphere  for  a  given  altitude  (cm). 

Subroutine  sckt,  adapted  from  LAMB,  is  shown  in  Figure  A-3.  It  computes  scale  factors  for  time, 
distance  and  pressure  for  use  with  the  1-KT  Standard  given  the  altitude  (cm)  and  yield  (kt)  of  the  burst 
and  the  altitude  (cm)  of  the  station  of  interest.  Unneeded  scale  factors  have  been  eliminated  from  this 
routine. 

Function  dist  is  shown  in  Figure  A-4.  It  computes  the  radius  from  a  stack  to  a  station  given  their 
coordinates. 

Function  tshk  is  shown  in  Figure  A-4.  It  computes  the  time  of  arrival  of  the  shock  wave  produced 
by  a  1-KT  burst  given  the  scale  radius  from  the  center  of  burst  and  scales  it  to  the  yield  of  the  specified 
burst. 

Function  dtpp  is  shown  in  Figure  A-4.  It  computes  the  positive  phase  duration  of  the  blast  wave 
produced  by  a  1-KT  burst  given  the  scale  radius  from  the  center  of  burst  and  scales  it  to  the  yield  of  the 
specified  burst. 

Function  peak  is  shown  in  Figure  A-5.  It  computes  the  peak  overpressure  of  the  blast  wave  associated 
with  a  specified  burst  given  the  scale  radius  from  the  center  of  burst. 
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subroutine  at62 (tty, wsp, cs, wsr, wst) 

computes  standard  atmosphere  conditions  at  a  given  altitude 
(adapted  from  DNA  1-KT  STANDARD) 

tabat(l)  =  r,  the  gas  constant  in  erg/mole/K 
tabat(2)  =  radius  of  the  earth  in  cm  ,  .  ,  , 

tabatO)  =  acceleration  due  to  gravity  at  sea  level  in  cm/s/s 
tabat(4)  =  molecular  weight  of  air  at  sea  level 


tabz  =  altitude  in  cm 

tabt  =  molecular  scale  temperature  in  K, 

tabl  =  molecular  scale  temperature  gradient  in  K/cm 

tabp  =  pressure  in  dyne/cm^*2 

nz  =  the  number  of  altitudes 

dimension  tabat (4 ) , tabz (22) , tabl (21 ) , tabt (22) , tabp (22 ) 

data  for  temperate  atmosphere  _ , 

data  nz/  21/,  rhoz/  1 .22500000e-03/, 

+  8.3144000e+07, 6.367488e+08,9.80665000e+02,2.8964 
data  tabz/ 

0.  ,  1.10190000e+06, 

-  ■■  5.24290000e+06, 

l.OOOOOOOOe+07, 

1.60000000e+07, 

3.00000000e+07, 

7.00000000e+07/ 


4 .73500000e+06 
9.00000000e+06 
1.50000000e+07 
2.20000000e+07 
6.00000000e+07 
data  tabl/ 
+-6.49291769e-05 
+-l,72248474e-07 
+  2.98166740e-05 
+  1.49589024e-04 
+  3.490350006-05, 
+  1.091624206-05/ 
data  tabt/ 

+  2.88150000e+02 
+  2.707041376+02 
+  1.80736048e+02 
+  9.60857386e+02 
+  1.551014046+03 
+  2.59146286e+03 
data  tabp/ 

+  1.01325D00e+06 
+  1.10904998e+03 
+  1.64379881e+00 
+  5.06169601e-03 
+»8.67381898e-04 
+  3.55894454e-06 


2.00630000e+06, 

6.15910000e+06, 

l.lOOOOOOOe+07, 

1.700000006+07, 

4.000000006+07, 


tabat/ 

400e+01/ 

3.21620000e+06, 

7.99940000e+06, 

1.20000000e+07, 

1.90000000e+07, 

5.000000006+07, 


9.28018576e-08,  9. 86255062e-06, 
-1. 96000240e-05,-3.91696897e-05, 
5.020200806-05,  9 . 97762300e-05 , 


6.97598500e-05, 

2.58868500e-05, 


1.004074906-04, 

3.31099360e-05, 

2. 16604540e+02,  2.16688470e+02, 
2.70616652e+02,  2.52659110e+02, 
2.10552722e+02,  2.60754730e+02, 
1.11044641e+03,  1 .21085390e+03, 
1.83024204e+03,  2. 16134140e+03, 
2.700625286+03/ 

2.26320000e+05,  5 . 47486994e+04 , 
5.90004987e+02,  1.82098959e+02, 
3.007497816-01,  7.35439451e-02, 
3.69429709e-03,  2.79259780e-03, 
1.943174306-04,  4 . 15743157e-05, 
1.22936355e-06/ 


2.77080373e-05, 

1.60821507e-07, 

2.001088096-04, 

6.68801467e-05, 

1.71252960e-05, 


2.28621170e+02, 

1.805751306+02, 

3.60530960e+02, 

1.35037360e+03, 

2.420209906+03, 


8.68013979e+03, 

1.03769924e+01, 

2.52169805e-02, 

1.68519867e-03, 

1.13023464e-05, 


data  cons/l .389779604el4/ 

find  the  table  index  corresponding  to  the  given  altitude 
if (tty . gt . 0 . 0)  then 
do  llO  jat=l,nz 

if (tty-tabz ( jat) )  120,130,110 
110  continue 
jat=nz+l 
120  1at=jat-l 

3at=max0 (1, jat) 
else 
jat=l 
endif 

compute  the  standard  atmosphere  conditions 
130  dum2= (tabz ( jat) -tty) / (tabat (2) +tty) / (tabat (2) +tabz ( jat) ) 
dum3=(tabat (2)+tabz  ^at) ) / jtabat (2) +ttY) 
varl= 
var2= 
f s=-c _ 

wsp=tabp ( jat) *exp (f s)  . 

wst=tabt ( jat ) +tabl ( jat ) * ( tty-tabz ( ^at ) ) 
wsr=wsp*tabt (1) *rhoz/ (wst*tabp(l) ) 
cs=sqrt (1 .4*wsp/wsr) 

return 

end 
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Figure  A-2.  Subroutine  at62  from  the  waveform  module. 
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c _ _ _ _ 

c 

subroutine  sckt (zbur, wbur, zsta,  tscal, rscal,pscal) 
c 

c  confutes  scale  factors  for  use  with  1-kt  standard 

C  (adapted  from  DNA  1-KT  STANDARD) 

c  independent  variables : 

c  zbur  -  burst  altitude  in  cm  above  sea  level 

c  zsta  *  station  altitude  in  cm  above  sea  level 

c  wbur  =  yield  of  the  burst  in  kt 

c 

c  dependent  variables : 

c  tscal  ■=  scales  actual  times  to  1-kt  sea  level 

c  rscal  =  scales  1  kt  sea  level  dimensions  to  actual 

c  pscal  =  scales  1  kt  sea  level  pressures  to  actual 

c 

common/conf/  thrd, tpi 
common/cona/  pl,cl,rl,tl 
c 

call  at62  (zsta,  p3,c3,r3,t3) 
pscal=p3/pl 

if ( (zsta-zbur) .ge . 100 . )  then 
call  at62(zbur,  p2,c2,r2,t2) 
dz=zsta-zbur 
b=alog (p3/p2) /dz 
bzsta=b*zsta 
a=p3*exp (-bzsta) 

p3=a* (exp (bzsta) -exp {b*zbur) ) / (b*dz) 
b=alog(r3/r2) /dz 
bzsta=b*zsta 
a=r3*exp (-bzsta) 

r3=a* (exp (bzsta) -exp (b*zbur) ) / (b*dz) 
c3=sqrt (1 . 4*p3/r3) 
endif 

rscal= (wbur*pl/p3) **thrd 
tscal=c3/ (rscal*cl) 


end 
c  _ 
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Figure  A-3.  Subroutine  sckt  from  the  waveform  module. 
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c 

c 

c 


function  dist (xstan,ystan, zstan,xstck,ystck, zstck) 

computes  the  radius  from  a  stack  to  a  station 

independent  variables: 
xstan  =  x-coord^nate  of  the  station 
ystan  =  y-coordinate  of  the  station 
zstan  =  z-coordinate  of  the  station 
xstck  =  x-coordinate  of  the  stack 
ystck  =  y-coordinate  of  the  stack 
zstck  =  z-coordinate  of  the  stack 

^deftx'^^  x-comp6nent  of  separation  between  station  and  stack 

deity  =  y-component  of  separation  between  stat4.on  and  stack 

deltz  =  z-component  of  separation  between  station  and  stack 

deltx=xstan-xstck 

delty=ystan-ystck 

dfst=sqrttdeltx*deltx+delty*delty+deltz*deltz) 

return 

end 


function  tshk (tscal, rslkt) 

computes  the  shock  arrival  time ^f or  a 
specified  burst  as  a  function  of  radius 

independent  variables: 
tscal  =  time  scale  factor 
rslkt  =  1  kiloton  scale  radius 

dependent  variably:  , 
tshk  =  shock  arrival  time 


local  variable:  ,  .  , 

tslkt  =  1  kiloton  scale  arrival  time 


common/cons/  r0s,c0s,d0s 

compute  the  1-kt  shock  arrival  time 
tslkt= (rslkt-dOs*tanh (rslkt/rOs) ) 7c0s 

scale  the  arrival  time  to  actual 
tshk=tslkt/ tscal 


return 

end 


function  dtpp (tscal, rslkt) 

computes  the  positive  phase  duration  for 
a  specified  burst  as  a  function  of  radius 

independent  variables: 
tscal  =  time  scale  factor  ,, 
rslkt  =  1-kiloton  scale  radius 

dependent  variable: 
dtpp  =  positive  phase  duration 


local  variable:  ,  ,  .  ,  o 

dtlkt  =  1-kiloton  scale  positive  phase  duration 

common/conp/  a0p,b0p2, cOp, alp,blp,rOp 

apply  the  lower  limit  to  the  radius 
rl=amaxl (rslkt, rOp) 

dtlk?=cOpMl.i-aip?rl**blp)^sqft(a!fo|uOp*sqrt  (rl*rl+b0p2) ) ) 


scale  the  positive  phase  duratioh  to  actual 
dtpp=dtlkt/tscal 

return 

end 
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Figure  A-4.  Functions  dist,  tshk,  and  dmp  from  the  wavefomi  module. 
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function  peak (pscal, rsllb) 

computes  the  peak  overpressure  as  a  function 
of  scaled  radius  using  the  kingery  fit 

independent  variables : 
pscal  =  pressure  scale  factor 
rsllb  =  1-pound  scale  radius 

dependent  variable: 
peak  «=  peak  overpressure 

common/conv/  okt, ocm, odsc, ocmkt3 
common /conk/  ak0,ak(8) 

compute  the  peak  overpressure  in  psi 
alr^=alog  (rsllb) 
ping=akO 
do  110  i=l,8 

ping=ping+ak (i) *alrsl**i 
110  continue 

ping=exp (ping) 

convert  to  dyne/cm**2  and  scale  to  sea  level 
peak=odsc*pscal*ping 

return 

end 


function  reak (prtio, rambi) 

corrputes  peak  density  (adapted  from  LAMB) 

independent  variables; 
prtio  =  pressure  ratio 
rambi  =  ambient  density 

dependent  variable 
reak  =  peak  density: 

common/coni/  gamma, gmone,gpone 

reak=rambi* ( (gpone*prtio+gmone) / (gmone*prtio+gpone) ) 

return 

end 


c 

c 

c 


function  veak (ppeak, rpeak, rambi) 

computes  the  peak  particle  velocity  (adapted  from  LAMB) 

independent  variables : 
rambi  =  ambient  density 
ppeak  =  peak  overpressure 
rpeak  =  peak  overdensity 

dependent  variable:  ,  x 

veak  •=  peak  particle  velocity 

opeak=rpeak- rambi 

veak=sqrt (ppeak*opeak/ (rambi* (rambi+opeak) ) ) 

return 

end 
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Figure  A-5.  Functions  peak,  reak,  and  veak  from  the  waveform  module. 
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Function  reak,  adapted  from  LAMB,  is  shown  in  Figure  A-5.  It  computes  the  peak  density  of  the 
Wast  wave  associated  with  a  specified  burst  given  the  pressure  ratio  and  the  ambient  density.  It  is  only 
used  in  conjunction  with  the  LAMB  peak  overpressure  combination  algorithm. 

Function  veak,  adapted  from  LAMB,  is  shown  in  Figure  A-5.  It  computes  the  peak  particle  velocity 
of  the  blast  wave  associated  with  a  specified  burst  given  the  peak  overpressure  and  overdensity,  and  the 
ambient  density.  It  is  only  used  in  conjunction  with  the  LAMB  peak  overpressure  combination  algorithm. 

COALESCENCE  MODULE 

These  routines  determine  the  effective  yield  at  each  station  for  a  specified  initial  donor  stack.  It  is 
recognized  that  more  than  one  group  of  coalesced  waves  may  pass  a  station  and  the  yield  selected  is  that 
associated  with  the  group  that  produces  the  largest  combined  peak  overpressure. 

Subroutine  wave  is  shown  in  Figure  A-6.  It  cycles  through  the  stacks,  computes  the  required  scale 
factors  with  calls  to  sckt,  computes  the  stack-to-station  radii  with  calls  to  dist,  and  computes  the  waveform 
variables  with  calls  to  tshk,  dtpp,  and  peak,  as  well  as  to  reak  and  \eak  when  the  LAMB  model  is 
specified. 

Subroutine  ssup  is  shown  in  Figure  A-7.  It  totals  the  yields  and,  using  simple  superposition,  combines 
the  peak  oveipiessuies  associated  with  coalesced  blast  waves.  This  option  can  only  be  selected  by  altering 
the  default  selection  in  the  file  bwaco.cfg. 

Subroutine  lamb  is  shown  in  Figure  A-7.  It  totals  the  yields  and,  using  the  LAMB  algorithm, 
combines  the  peak  overpressures  associated  with  coalesced  blast  waves.  This  option  can  only  be  selected 
by  altering  the  default  selection  in  the  file  bwaco.cfg. 


Subroutine  comb  is  shown  in  Figure  A-7.  It  totals  the  yields  and,  usingjhe  total  jield,  computes  the 
peak  overpressures  associated  with  coalesced  blast  waves. 

Subroutine  coal  is  shown  in  Figure  A-8.  It  is  given  the  index  of  a  specified  stack  and  the  coordinates 
of  the  current  station.  It  cycles  through  all  the  other  stacks  and,  using  the  shock  and  zero-overpressure 
arrival  times  for  all  stacks  at  the  current  station,  it  applies  the  coalescence  criterion  that  the  shock  from 
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subroutine  wave (istan, jstan) 

cycles  through  the  stacks  and  computes  the  waveform  variables 

independent  variables: 
istan  =  station  x-^ndex 
jstan  =  station  y-index 

configuration  variable;  ,  ,  .  ^  . 

icomB  =  pressure  combination  algorithm  switch 

station  variables:  ^  ^ 

xstat  =  array  of  x-coordinates  for  the  stations 
ystat  =  array  of  y-coordlnates  for  the  stations 
zstat  =  z-coordinate  of  stations 


stack  variables: 
nstak  =  total  number  of  stacks 
xstak  =  array  of  x-coordf nates  for  the  stacks 
ystak  =  array  of  y-coordinates  for  the  stacks 
zstak  =  z-coordlnate  for  the  stacks 
wstak  =  array  of  yields  for  the  stacks 
tstak  =  array  of  stack  initiation  times 

tshok  =  array  of  shock  arrival  times  for  wav^s  from  stacks 
tzero  =  array  of  zero-overpressure  arrival  times 
for  waves  from  stacks 

ppeak  =  array  of  peak  overpressures  for  waves  from  stacks 
rpeak  =  array  of  peak  densities  for  waves  from  stacks 
vpeak  =  array  of  peak  velocities  for  waves  from  stacks 

local  variables: 
istak  =  stack  counter 
tscal  =  time  scale  factor 
rscal  =  radius  scale  factor 
pscal  =  pressure  scale  factor 
rstst  =  stack  to  station  radius 
rslkt  =  radius  scaled  to  1  kiloton 
rsllb  =  radius  scaled  to  1  pound 
prtio  =  pressure  ratio 


common/conv/ 
common/icfg/ 
common/ambi/ 
common/sttp/ 
common/stkn/ 
common/stkp/ 
common/ St kt/ 
common/stky/ 
common/stkw/ 


okt, ocm, odsc, ocmkt3 
icomb, iscrn, iprnt 
pambi, rambi 

xstat  (120) , ystat (120) , zstat 
nstak 
xstak 
vstak 


wstak 

tshok 

ppeak 


ioo 


100 

100 

100 


,  tzero (100) , 

, rpeak (100) , vpeak (100) 


c 

c 

c 

c 


cycle  through  the  stacks 
do  110  istak=l, nstak 

compute  the  scale  factors 

call  sckt (zstak, wstak  (istak) , zstat,  tscal, rscal, pscal) 

compute  the  stack  to  station  radius 
rstst=dist (xstat (istan) , ystat (jstan) , zstat, 
xstak (istak) , ystak (istak) , zstak) 

scale  the  radius  to  1-kt 
rslkt=rstst/rscal 

compute  the  shock  and  zero-overpressure  arrival  times 
tshok (istak) =tshk (tscal, rslkt ) +tstak (istak) 
tzero (istak) =dtpp (tscal, rslkt) +tshok (istak) 

scale  the  radius  to  1-lb  and  nondoubled  yield 
rslIb=ocmkt3*rslkt 


compute  the  peak  overpressure,  density  and  particle  velocity 
ppeak (istak )=peak (pscal, rsllb) 
if  (icomb, eg.  ll  then, _ _ _  , _ 


prtio= ippeax (5stak) +p^mbi) /p^n*! 
rpeak (istak)=reak (prtip,fambi) 
vpeak (Istak) =veak (ppeak (Istak) , rpeak (istak) , rambi) 
rpeak (Istak) =rpeak (Istak) -rambi 


rambi) 


endif 


110  continue 


return 

end 
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Figure  A-6.  Subroutine  wave  from  the  coalescence  module. 
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subroutine  saup (ncoal, wcoal, pcoal^  wconb, pconb) 

COTbinea  yields  and  peak  overpreaaurea  using  ainple  superposition 


independent  ya 


wcoa 

pcoa 


i: 


riables 
__r  of  CO. 

SI 


elds  of  sources  of  coalesced  waves 
ak  overpressures  fron  coalesced  waves 


dependent  var  .ables: 
wcoi^  •  combined  yiel< 
omb  ■>  confined  ^ak 


pcomb  ■>  combined  peak  overpressure  o: 

local  variable: 
iwave  ••  wave  counter 

dimension  wcoal (100) ,pcoal (100) 


d  from  sources  of  all  coalesced  waves 
-)f  all  coalesced  waves 


wcomb-Q. 

8  comb* 0 , 

0  110  iwave«l, ncoal, 
wcomb-wcomMwcoal  (4-vave) 
pcomb-pcomb-i-pcoal  uwave) 
110  cofitinue*^ 

return 

end 


subroutine  lasib  (ncoal ,  rasta, wcoal « pcoal ,  rcoal ,  vcoal^  wcomb»  pconb) 
e^Bbines  yields  and  peak  overpressures  using  IAMB  algorithm 
independent  variables: 

ncoal  -  nui^r  of  coalesced^waves  ^  ^  . 

rasta  *  ambient  density  at  the  current  station 

'  elds  of  sources  o^  coalesced  waves 


pcoai  •  arra9  of  beak  overpressures  from  coalesced  waves 
rcoal  *  array  of  Mak  densities  from  coalesced  waves 
vcoal  •  array  oi  ^ak  particle  velocities  from  coalesced  waves 

dependent  variables:  > 

wcomb  -  combined  yield  from  sources  of  all  coalesced  waves 
pcomb  *  cOTbined  peak  overpressure  of  all 


SI 


seed  wav 
alesced 


pcomb 

local  variables: 
iwave  *  wave, counter 


coalesced  waves 


waves 

_  waves 

coalesced  waves 


dimension  wcoal (100) , pcoal (100) « rcoal (100) , vcoal (100) 

wcomb*0 . 
rcomb*rasta 
pcomb-0 . 
vcomb-Q. 
ecomb*0 , 

do  110  iwave-1, ncoal 

weomb-wcomb-f wcoal  (iwave) 
rcomb*rcomb+ rcoal (iwave)  ,  . 
rv-(rasta'^rcoal  (iwave) )  ^vcoal  (iwave) 
vcomb-vcomb+rv  .  , 
ecomb*ecomb‘t'rv*vcosl  (iwave) 
pcomb-pcombtpcoal (iwave) 

110  continue 

^om^^orab+5?S*  (ecomb-rcomb*vcoinb*vcoffib) 

return 

end 

subroutine  comb ( ncoal, xcoal,ycoal, wcoal, X8tat,ystat, sstat, wcomb, 
*  pcomb) 


cl  *  scale  factor  for  combined  stacks 

cl  •  radius  scale  factor  for  combined  stacks 
cl  *  pressure  scale  factor  for  combined  stac. 


_  stacks 

combined  stacks 
combined  stacks 


pcoi 


ite  the  combined 
^ak(pcscl,rsl* 


y* 

ressure 


zstak) /reseX 
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combines  peak  overpressures  using  pressure  from  combined  stacks  ^§ooS 

independent  variables:, 
ncoal  -  number  of  coalesced  waves 

xcoai  -  array  of  x-coordinates  of  sources  of  coalesced  waves 

«coaI  *  array  of  y^coordinates  of  sources  of  coalesced  waves 
coal  •  array  of  yields  of  sources  of  coalesced  waves 
xstat  *  X'coordinate  of  the  current  station 
ystat  -  y-coordinate  of  the  current  station 
tstat  -  t-coordinate  of  the  current  station 

dependent  variables :  .  , ,  ,  ^ 

wcomb  *  combined  yield  from  sources  of  all  coalesced  waves 
pcomb  ••  combined  ^ak  overpressure  of  all  coalesced  waves 

stack  variable:  ,  ,  ^ 

zstak  -  z-coordinate  for  the  stacks 

local  variables: 
tcsci 
rcscl 

pcsci  -  pressure _ _ 

xcomb  *  x-coordinate  of  center 
ycomb  •  y-coordinate  of  center 
rsllb  -  1-lb  scale  radius 

conmon/conv/  okt,ocm.odsc,pcmkt3,,,  . 
comreon/stkp/  xstak (100) , ystak(lOO) , zstak 

dimension  xcoai (100) , ycoal <100) , wcoal (100) 

find  weight  and  center  of  combined  stacks  •  . 

call  sseqlncoal, xcoai, ycoal, wcoal,  xcomb, ycomb, wcomb) 

compute  the  scale  factors  ^ 
call  sekt (zstak. wcomb, zstat,  tcsci, rcscl, pcsci) 


g888i 


return 

end 


Figure  A-7.  Subroutines  ssup,  lamb,  and  comb  from  the  coalescence  module. 
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subroutine  coel <lepec»xetat.yetet» xstet, rests,  ncoel, wcomb,pcoiib} 

applies  the  coelescence  criterion  end  detemlnes  the 
tetel  yield  end  conblned  peak  overpressure  from  ell 
neves  which  coalesce  with  the  weve  from  e  specified  stack 

independent  variables: 
ispec  -  index  of  the  specified  stack 
xstat  -  x-coordinate  oi  the  current  station 
ystat  -  y-coordinate  of  the  current  station 
kstat  *  {-coordinate  of  the  current  station 
rasta  *  ambient  density  at  the  current  station 


dependent  variables 
ncoal  •  nui^r  of 


_  _ ,  _  -  _  ,  lalesced  waves 

wcomb  •  ccmiblned  yield  of  the  sources  of  all  eolaesced  waves 
pc«&b  -  combined  peak  overpressure  of  all  coalesced  waves 


configuration  variable:  ,  ^  ^ 

ieomb  »  pressure  combination  algorithm  switch 

stack  variables: 
nstak  -  total  number  of  stacks 
xstak  •  array  of  x-coordlnetes  for  the  stacks 
ystak  •  array  of  y-coordinates  for  the  stacks 
wstak  »  array  of  yields  for  the  stacks 

ppeak  •  array  of  peak  overpressures  for  waves  from  the  stacks 
rnak  array  of  Msk  densities  for  waves  from  the  stacks 
VMak  -  array  of  peak  velocities  for  waves  from  the  stacks  ^ 
tshok  •  array  of  shock  arrival  times  for  waves  from  the  stacks 
tsero  •  array  of  xero-overpressure  arrival  times 
for  waves  from  the  stacks 

local  variables: 
istak  •  stack  counter 

xcoal  •  array  of  x-coordinates  of  sources  of  coalesced  waves 

ycoal  •  array  of  y-coordlnates  of  sources  of  coalesced  waves 

wcoal  •  array  of  yields  of  sources  of  coalesced  waves 

pcoal  *  array  of  peak  overpressures  from  coalesced  waves 

rcoal  •  array  of  ^ak  densities  from  coalesced  waves 

vcoal  •  array  of  peak  particle  velocities  from  coalesced  waves 

common/lcfq/  ieomb, Iscrn, Iprnt 
common/stkn/  nstak 

common/ St kp/  xstak (100) , ystak ( 100) , xstak 

comnon/stky/  wstak(lOO) 

common/ St kw/  tshok (100) ,t zero (100) , 

4  ppeak I 100) , rpeak UOO) , vpeak(lOO) 

dimension  xcoal (100) , ycoal (100) , wcoal (100) , 

4  l^oal  (100) ,  rcoal  (100) ,  vcoal  (100) 

Initialize  the  number  of  coalesced  waves  and  coalescence  arrays 
ncoal-1 

xcoal (1) -xstak (ispec) 
ycoal (1 I -ystak (ispec) 
wcoal ( 1 ) -wstak (ispec) 
pcoal ( 1 ) -ppeak (ispec) 
if (ieomb. eq. 1 )  then 
rcoal (l)-rpeak (ispec) 
v^oal ( 1 ) -vpeak (ispec) 

cycle  through  the  stacks 
do  110  istak-1, nstak 

if  the  current  stack  is  not  the  specified  stack 
if (istak. ne. ispec)  then 

if  the  wave  from  the  current  stack  is  coalesced 
with  the  wave  from  the  specified  stack 
if  jtshok (istak)  .^e.tshokdspec) )  then 


if (tshok( 


. ie.tzero (Ispec) )  then 


increment  the  number  of  coalesced  waves 
ncoal-ncoal4i 

store  the  peak  overpressure,  overdensity  and  velocity 
xcoal (ncoal) -xstak (istak) 
ycoal (ncoal ) -ystak (xstak) 
wcoal (ncoal ) -wstak (istak) 
pcoal (ncoal) -ppeak (istak) 
if (ieomb. eq.l)  then 

rcoal (ncoal) -rpeak (istak) 
v^oal (ncoal) -vpeak (istak) 


endif 

endlf 


endif 

110  continue 


combine  the  yields  and  peak 
overpressures^  of  the  coalesced  waves 


if (ncoal. gt. 1)  then 

use  sliule  superposition 

if (ieomb. eq.O)  call  ssup(ncoal, wcoal, pcoal,  wcomb, pcomb) 

use  the  IAMB  algorithm  ,  «  .  , 

if (ieomb. eq. 1)  call  lamb(ncoal,r«sta, wcoal, pcoal, rcoal, vcoal. 


wcoiBb.  pcomb) 

use  the  combined  dtacks 

if (ieomb. eq. 2)  call  bomb(ncoal,xcoaI. ycoal, wcoal* xstat, ystat, 
4  zstat,  wcomb, pcomb) 

else 

use  the  yield  and  peak  overpressure  from  the  specified  stack 
wcomb-wcoal (1) 
peomb-pcoal(l) 

endif 

return 

end 
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Figure  A-8.  Subroutine  coal  from  the  coalescence  module. 
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wie  wave  lie  in  tlw  positive  overpressure  phase  of  another  wave.  In  this  way,  it  determines  the  number 
of  stacks  pnxiucing  waves  which  coalesce  with  the  wave  from  the  spedfied  stack  as  well  as  the  total  yield 
aiKl  combined  peak  overpressure  (using  calls  to  ssup,  lamb,  or  comb  depending  on  how  the  session  is 
configured)  of  the  coalesced  waves.  It  returns  the  number  of  coalesced  waves  as  well  as  the  associated 
total  yield  and  combined  peak  overpressure. 

Subroutine  is  shown  in  Figure  A-9.  It  is  given  the  coordinates  of  flie  current  station.  It  cycles 
through  the  stacks,  calling  coal,  and  determines  the  effective  yield  associated  with  the  group  of  coalesced 
waves  producing  the  largest  combined  peak  oveipressure  at  the  current  station.  It  returns  the  number  of 
coalesced  waves  in  the  group,  the  effective  yield  and  the  combined  peak  overpressure. 

Subroutine  ^  is  shown  in  Figure  A- 10.  It  cycles  through  the  stations  and  computes  the  waveform 
variables  for  each  stack  with  calls  to  wave  and  the  effective  yields  and  combined  peak  overpressures  with 
calls  to  y^.  It  saves  the  peak  overpressure  and  effective  yield  associated  with  the  largest  combined  peak 
overpressure  at  each  statioa 

PROPAGATION  MODULE 

These  routines  detennine  the  order  in  which  the  stacks  detonate  and  the  time  of  detonation  of  each 
stack  for  a  given  initial  donor. 

Function  tdly  is  shown  in  Figure  A-1 1.  It  is  given  the  coordinates  of  a  donor  and  an  acceptor  stack 
as  well  as  the  propagation  velocity  of  the  donor.  It  computes  and  returns  the  initiation  delay 
(communication  time)  between  the  donor  and  acceptor  stack. 

Subroutine  arty  is  shown  in  Figure  A- 12.  It  is  given  the  coordinates,  initiation  time,  and  propagation 
velocity  of  a  donor  stack.  It  cycles  through  all  the  other  stacks  which  have  not  yet  become  donors  and, 
with  calls  to  tdly,  determines  flie  initiation  time  of  the  earliest  acceptor  for  tiiat  donor.  It  assigns  that  time 
to  the  stack  initiation  time  array  if  it  is  earlier  than  the  time  already  assigned.  It  returns  the  index  of  that 
acceptor  and  its  initiation  time. 
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subroutine  yef f (xstat, ystat, zstat, rambi,  icsta, wstat,pstat) 

computes  the  effective  yield  and  combined 
peak  overpressure  at  a  current  station 

independent  variables : 

xstat  =  x-coordinate  of  the  current  station 
ystat  =  y-coordinate  of  the  current  station 
zstat  =  z-coordinate  of  the  current  station 
rambi  =  ambient  density  at  the  current  station 

dependent  variables; 

icsta  =  coalescence  status  at  the  current  station 

wstat  =  effective  yield  at  the  current  station 

pstat  =  combined  peak  overpressure  at  the  current  station 

stack  variable: 

nstak  =  total  number  of  stacks 

local  variables: 
istak  =  stack  counter 

ncoal  =  number  of  waves  coalesced  with 
the  wave  from  the  current  stack 
wcomb  =  combined  yield  at  the  current  station 
pcomb  =  combined  peak  overpressure  at  the  current  station 

common/stkn/  nstak 

initialize  the  coalescence  status,  effective 
yield  and  combined  peak  overpressure 
icsta=0 
wstat^O . 
pstat=0 . 

« 

cycle  through  the  stacks 
do  110  istak=l, nstak 

determine  the  total  yield  and  combined  peak  overpressure  from 
all  waves  which  coalesce  with  the  wave  from  the  current  stack 
call  coal (istak, xstat, ystat, zstat, rambi,  ncoal, wcomb, pcomb) 

save  the  coalescence  status,  total  yield 
and  combined  peak  overpressure  associated 
with  the  largest  combined  peak  overpressure 
if (pcomb. gt .pstat)  then 
icsta=ncoal 
wstat=wcomb 
pstat=pcomb 
endif 

110  continue 
return 
end 
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Figure  A-9.  Subroutine  yeff  from  the  coalescence  module. 


61 


oo  o  o  noon  non  on  no  on  on  on  o o o o o o o o o o o o o o o o o o o o o o o o  on 


subroutine  effy 

computes  the  effective  yield  associated  with  the 
largest  combined  peak  overpressure  at  each  station 

station  variables: 
nstat  =  number  of  stations 
Iblok  =  length  of  station  blocks 
nblok  =  number  of  station  blocks 
xstat  =  array  of  x-coordinates  for  the  stations 

ystat  =  array  of  y-coordinates  for  the  stations 

zstat  =  z-coordinate  of  stations 

icoal  =  array  of  coalescence  statuses  for  the  stations 

wstat  =  array  of  effective  yields  for  the  stations 

pstat  =  array  of  effective  peak  overpressures  for  the  stations 
rambi  =  ambient  density 


local  variables: 
istat  =  station  counter 
istan  =  station  x-index 
jstan  =  station  y-index 

iccur  =  coalescence  status  at  current  station 

wscur  =  effective  yield  at  current  station 

pscur  =  combined  peak  overpressure  at  current  station 


common/alow/  pallo 
common/ambi/  pambi, rambi 
common/sttn/  nstat 
common/sttb/  Iblok, nblok 

common/sttp/  xstat ( 120) , ystat (120) .zstat 
common/ sttv/  icoal (3600) ,pstat(3600) , wstat (3600) 


initialize  the  station  report 
write  (6, 2010)  0, nstat 

cycle  through  the  stations 
do  110  istat=l, nstat 

report  the  current  station 
write  (6, 2020)  istat, nstat 

compute  the  station  indices 
istan=mod (istat-1 , Iblok) +1 
jstan= (istat-1 ) /lblok+1 

compute  waveform  variables 
call  wave (istan, jstan) 

compute  the  effective  yield  and  combined 
peak  overpressure  at  the  current  station 
call  yeff (xstat (istan) , ystat ( jstan) , zstat, rambi, 
+  iccur, wscur, pscur) 

save  the  effective  yield  and  combined  peak 
overpressure  associated  with  the  largest 
combined  peak  overpressure  for  any  donor 
if (pscur . gt .pstat  (istat) )  then 
if (pscur . It .pallo)  then 
icoal (istat ) =0 
wstat  (istat ) =0 . 

©ISG 

icoal (istat) =iccur 
wstat (istat)=wscur 
endif 

pstat (istat ) =pscur 
endif 

110  continue 


return 


2010 

2020 

c 


c 


format  statements 

f ormat ( '  Computing  coalescence  at  station 

format ('+  Computing  coalescence  at  station  ' 

end 


of 

of 


M4) 

'  ,14) 
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Figure  A- 10.  Subroutine  efiy  from  the  coalescence  module. 
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function  tdly (xacpt , yacpt , xdonr, ydonr, vdonr) 

computes  the  initiation  delay  between  a  given  donor  and  acceptor 

independent  variables : 
xacpt  =  x-coordinate  of  the  acceptor 
yacpt  =  y-coordinate  of  the  acceptor 
xdonr  =  x-coordinate  of  the  donor 
ydonr  =  y-coordinate  of  the  donor 
vdonr  =  propagation  velocity  of  the  donor 

dependent  variable: 
tdly  =  initiation  delay 

local  variables: 
deltx  =  x-coordinate  difference 
deity  =  y-coordinate  difference 

if (vdonr .eq. 0 . 0)  then 

infinite  propagation  velocity 
tdly=0 . 0 

else 

finite  propagation  velocity 
deltx=xacpt-xdonr 
delty=yacpt-ydonr 

tdly=sqrt (deltx*deltx+delty*delty) /vdonr 

endif 

return 

end 
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Figure  A-11.  Function  tdly  from  the  propagation  module. 
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subroutine  arly (xdonr, ydonr, tdonr, vdonr,  ierly,terly) 

computes  earliest  acceptor  for  a  given  donor 

independent  variables : 
xdonr  =  x-coordinate  of  donor 
ydonr  =  y-coordinate  of  donor 
tdonr  =  initiation  time  of  donor 
vdonr  =  propagation  velocity  of  donor 

stack  variables: 

xstak  =  array  of  x-coordinates  for  the  stacks 
ystak  =  array  of  y-coordinates  for  the  stacks 
tstak  =  array  of  initiation  times  for  the  stacks 
mstak  =  array  of  initiation  statuses  for  the  stacks 

dependent  variables: 
ierly  =  index  of  earliest  acceptor 
terly  =  initiation  time  of  earliest  acceptor 

local  variable: 
istak  =  stack  counter 

common/stkn/  nstak 

common/stkp/  xstak (100) , ystak (100) , zstak 
common/ stkt/  vstak  (100) , tstak (100) 
common/inst/  mstak (100) 

initialize  earliest  initiation  time 
terly=l.e9 

cycle  through  the  stacks 
do  110  istak=l, nstak 

if  the  current  stack  is  not  already  a  donor 
if (mstak (istak) .eq. 0)  then 

compute  its  initiation  time  with  respect  to  given  donor 
tstak (istak) =tdonr+tdly (xstak (istak) ,  ystak (istak) , xdonr, ydonr, 
+  vdonr) 

save  the  earliest  initiation  time 
if (tstak (istak) . It .terly)  then 
terly=tstak (istak) 
ierly=istak 
endif 

endif 

110  continue 
return 
end 
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Figure  A- 12.  Subroutine  arly  from  the  propagation  module. 
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Subroutine  prop  is  shown  in  Figure  A-13.  It  is  given  the  index  of  the  current  initial  donor  stack.  It 
cycles  through  all  the  stacks  which  have  become  donors  to  determine  the  next  donor  using  calls  to  arly. 
It  repeats  the  process  using  that  donor  and  continues  until  initiation  times  have  been  assigned  to  aU  of  the 
stacks. 

EXECUTION  MODULE 

These  routines  read  the  user’s  input,  set  up  the  requested  computation,  and  execute  it 

Subroutine  inpt,  shown  in  Figure  A-14,  reads  the  user’s  input  from  the  file  bwaco.inp. 

Subroutine  cent,  shown  in  Figure  A-15,  computes  either  coordinate  of  the  center  of  charge  of  a  given 
stack  arrangement. 

Subroutine  sseq,  shown  in  Figure  A-15,  calls  cent  to  compute  both  coordinates  of  the  center  of  charge 
of  a  given  stack  arrangement. 

Subroutine  stak,  shown  in  Figure  A-16,  calls  sseq  to  compute  the  single-stack  equivalent  of  all  the 
stacks,  converts  the  stack  variables  to  1-KT  Standard  units  and  doubles  their  yields  to  simulate  the  effects 
of  the  ground  plane. 

Subroutine  blok,  shown  in  Figure  A-17,  determines  the  station  arrangement  that  best  covers  the  region 
of  interest. 

Subroutine  stns,  shown  in  Figure  A- 18,  determines  the  region  of  interest  (if  none  is  specified  by  the 
user),  determines  the  station  arrangement  by  calling  blok,  defines  the  station  coordinates,  and  initializes 
the  station  variables. 

Subroutine  rtim,  shown  in  Figure  A- 19,  computes  and  reports  the  run  time  for  the  coalescence 
computation. 
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subroutine  prop(idoni) 

computes  initiation  times  of  stacks 

independent  variable: 

Idoni  =  index  of  initial  donor 

stack  variables: 
nstak  =  number  of  stacks 

xstak  =  array  of  x-coordinates  for  the  stacks 

ystak  =  array  of  y-coordinates  for  the  stacks 

vstak  =  array  of  stack  propagation  velocities 

tstak  =  array  of  initiation  times  for  the  stacks 
mstak  =  array  of  initiation  statuses  for  the  stacks 

local  variables: 
istak  =  stack  counter 

ierly  =  index  of  earliest  acceptor 

terly  =  initiation  time  of  earliest  acceptor 
idonx  =  index  of  next  donor 
tnext  =  earliest  acceptor  for  all  donors 
ndonr  =  number  of  identified  donors 

common/stkn/  nstak 

common/stkp/  xstak (100) , ystak (100) , zstak 
common/stkt/  vstak (100) , tstak (100) 
common/inst/  mstak (100) 

initialize  stack  initiation  statuses  and  times 
do  110  istak=l, nstak 
mstak  (lstak)=0 
tstak  (istak) =1 . e9 
110  continue 

initialize  donor  initiation  status  and  time 
ndonr=l 

mstak  (idoni) =1 
tstak (idoni ) =0 . 

find  earliest  acceptor  for  all  donors 
120  tnext=l.e9 

cycle  through  the  stacks 
do  130  istak=l, nstak 

if  the  current  stack  is  an  identified  donor 
if (mstak (istak) .gt . 0)  then 

determine  the  earliest  acceptor  for  that  donor 
call  arly (xstak (istak) , ystak (istak) , tstak (istak) , vstak (istak) , 
+  ierly, terly) 

save  the  earliest  acceptor  for  all  donors 
if (terly . It , tnext )  then 
idonx=ierly 
tnext =terly 
endlf 

endif 

130  continue 

Increment  donor 
ndonr=ndonr+l 

assign  earliest  acceptor  as  next  donor 
mstak (idonx) =ndonr 
tstak (idonx) =tnext 

check  for  completion  of  stack  initiation  time  assignments 
if (ndonr . It .nstak)  go  to  120 
c 
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Figure  A- 13.  Subroutine  prop  from  the  propagation  module. 
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subroutine  inpt 

obtains  problem  input  from  the  file  bwaco.inp 

problem  variable: 
ulabl  =  user' s  problem  label 

stack  variables: 
nstak  =  total  number  of  stacks 
xstak  =  array  of  x-coordinates  for  the  stacks 
ystak  =  array  of  y-coordinates  for  the  stacks 
zstak  =  z-coordinate  for  the  stacks 
wstak  =  array  of  yields  for  the  stacks 
vstak  =  array  of  stack  propagation  velocities 
idonr  =  array  of  initial  donor  switches 

region  variables: 

xmin  =  minimum  x-coordinate  of  region  of  interest 
xmax  =  maximum  x-coordinate  of  region  of  interest 
ymin  =  minimum  y-coordinate  of  region  of  interest 
ymax  =  maximum  y-coordinate  of  region  of  interest 

common/ulab/  ulabl 


1000010 
■  1000020 
1000030 
1000040 
1000050 
1000060 
1000070 
1000080 
1000090 
1000100 
1000110 
1000120 
1000130 
1000140 
1000150 
1000160 
1000170 
1000180 
1000190 
1000200 
1000210 
1000220 
1000230 
1000240 
1000250 
1000260 


xstak 

100) 

, ystak (100) , zstak 

1000270 

vstak 

100 

,tstak (100) 

1000280 

wstak 

100) 

1000290 

idonr 

100 

,  nidnr 

1000300 

character*80  ulabl, dummy 

log  call 
write(8,3010) 

open  the  user  input  file 

open (7, f ile=' bwaco . inp' , form=' formatted' , statu s= ' old', err=l 40) 

read  the  user' s  problem  label 
read (7,1010, end=lb0, err=150)  ulabl 
write(8,3020)  ulabl 

read  the  region  data 
read (7, 1010, end=l 60, err=150)  dummy 
read (7,1020,end=160, err=150)  xmin, xmax, ymin, ymax 
write (8, 3030)  xmin, xmax, ymin, ymax 

read  the  stack  data 
read (7, 1010, end=l 60, err=150)  dummy 
nstak=0 
zstak=0 . 

110  nstak=nstak+l 

if (nstak . gt . 100)  go  to  120 
read (7, 1020, end=120, err=150) 

+  xstak (nstak), ystak (nstak), wstak (nstak), vstak (nstak), idonr (nstak) 
write  (8,  3040) 

+  xstak (nstak) , ystak (nstak) , wstak (nstak) , vstak (nstak) ,  idonr (nstak) 
’  go  to  110 

close  the  user  input  file  and  adjust  the  number  of  stacks 
120  close(7) 

nstak=nstak-l 

remove  any  bogus  stacks 
130  continue 

if(xstak (nstak) .eq.O. . and .ystak (nstak) . eg. 0. . and . 

+  wstak  (nstak j . eg. 0 . . and. vstak (nstak) .eq.O. . and. 

+  idonr (nstak) .eg. 0)  then 
nstak=nstak-l 
go  to  130 
endif 

if (nstak. It . 1 )  go  to  170 
write (8, 3050)  nstak 
return 

errors  detected 

140  stop  '  BWACO:  ERROR  -  Unable  to  open  bwaco.inp!' 

150  stop  '  BWACO:  ERROR  -  Unable  to  read  bwaco.inp!' 

160  stop  '  BWACO:  ERROR  -  Premature  end  of  bwaco.inp!' 

170  stop  '  BWACO:  ERROR  -  Fewer  than  one  stack  specefied!' 

format  statements 
1010  format 
1020  format 
3010  format 


3020  format 
3030  format 
3040  format 
3050  format 

end 


a80) 
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inpt  called 
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Figure  A- 14.  Subroutine  inpt  from  the  execution  module. 
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subroutine  sseq (nstck, xstck, ystck, wstck,  xsseq, ysseq, wsseq) 

computes  the  single  stack  equivalent  of  the  distributed  stacks 

independent  variables: 
nstck  =  number  of  stacks 

xstck  =  array  of  x-coordinates  for  the  stacks 
ystck  =  array  of  y-coordinates  for  the  stacks 
wstck  =  array  of  yields  for  the  stacks 

dependent  variables : 

xsseq  =  x-coordinate  of  center  of  charge 
ysseq  =  y-coordinate  of  center  of  charge 
wsseq  =  single  stack  equivalent  weight 

local  variable: 
istck  =  stack  counter 

dimension  xstck (100) ,  ystck (100) , wstck (100) 

compute  the  single  stack  equivalent  weight 
wsseq=0 . 

do  110  istck=l, nstck 

wsseq=wsseq+wstck (istck) 

110  continue 

compute  the  center  of  charge 
xsseq=cent (nstck, xstck, wstck, wsseq) 
ysseq=cent (nstck, ystck, wstck, wsseq) 

return 

end 


function  cent (nstck, rstck, wstck, wtotl) 

computes  one  coordinate  of  the  center  of  charge 

independent  variables: 
nstck  =  number  of  stacks 

rstck  =  array  of  coordinates  for  the  stacks 
wstck  =  array  of  yields  for  the  stacks 
wtotl  =  total  yield  of  stacks 

local  variable: 
istck  =  stack  counter 

dimension  rstck (100) , wstck  (100) 

rwstk=0 . 

do  110  istck=l, nstck 

rwstk=rwstk+rstck (istck) *wstck (istck) 

110  continue 

cent=rwstk/ wtotl 

return 

end 
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Figure  A- 15.  Subroutine  ssea  and  function  cent  from  the  execution  module. 
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subroutine  stak 

defines  the  stacks 

stack  variables: 
nstak  *  total  number  of  stacks 
xstak  “  array  of  x-coordinates  for  the  stacks 
ystak  *  array  of  y-coordinates  for  the  stacks 
zstak  "  z-coordlnate  for  the  stacks 
wstak  =  array  of  yields  for  the  stacks 
vstak  •  array  of  stack  propagation  velocities 
nidnr  *  number  of  initial  donors  specified 

region  variables: 

xcent  *  x~coordinate  of  center  of  charge 

ycent  «=  y-coordinate  of  center  of  charge 

rinhb  *  inhabited  building  radius 

xmini  =  minimum  x-coordinate 

xmaxi  »  maximum  x-coordinate 

ymlnl  =  minimum  y-coordinate 

ymaxl  *  maximum  y-coordinate 

local  variables: 
istak  *■  stack  counter 
dregn  *  region  margin 
rregn  *  region  radius 

common/conf/  thrd,tpl 

common/ conv/  okt, ocm, odsc,ocmkt3 

common/regn/  xcent, ycent, rinhb, xmini, xmaxi, ymini,ymaxi 
common/alow/  pallo 
common/stkn/  nstak 

common/stkp/  xstak (100| , ystak (100) , zstak 
common/stky/  wstak (100) 
common/ St kt/  vstak (100) , t stak (100) 
common/stkd/  idonr  (100) , nidnr 
common/totl/  wtotl 

log  call 
write (8,3010) 

compute  single  charge  equivalent  and  inhabited  building  radius 
call  sseq  (nstak, xstak, ystak, wstak,  xcent, ycent, wtotl) 
if (wtotl?lt.3.0e4)  then 
if (wtotl. gt .100.)  then 
rinhb-380. 
else 

rinhb=204 . 
endif 

pallo=8 .3e4 
endif 

if (wtotl . ge .3 .0e4 .and.wtotl .It ,1. OeS)  then 
rinhb=l2, 2*wtotl**thrd 
pallo“8.3e4 
endif 

if (wtotl.ge.l .OeS.and.wtotl .lt.2.5e5)  then 
rinhb-0:738*wtotl**0.577 
pallo»9.7e4-0.14*wtotl 
endif 

if (wtotl.ge.2.5e5)  then 
rinhb=15.  H*wtotl**thrd 
pallo»6.2e4 
endif 

dregn-0 .05* rinhb 
rregn»rinhb+dregn 
xmini“X cent -rregn 
xmaxi-xcent+rregn 
ymini=ycent -rregn 
ymaxi-ycent+rregn 

write  (8,3020)  rinhb, xmini, xmaxi, yminl,ymaxi 

adjust  the  region  limits,  convert  the  stack  variables  to  1-kt 

standard  units  and  double  the  yields  for  ground  plane  simulation 
nidnr-0 
do  110  lstak=] 
rstak=15.24’ 
xmini=aminl 
xmaxi-amaxl 
ymini-aminl 
ymaxl -amaxl 
xstak (istak) 
ystak (istak 
wstak (istak) 
vstak  (istak) 
write  (8.3030) 

■K  xstak  (istak)  .ystak  (istak) , wstak  (istak) ,  vstak  (istak) 
if (idonr (istak) .ne.O)  nidnr«nidnr+l 
continue 


.stak) **thrd+dregn 
.stak) -rstak, xmini) 
.stak) +rstak, xmaxi) 
.stak j -rstak, ymini I 
.stak) +rstak,ymaxi) 
(istak) /ocm 

^ _ (istak)  /ocm 

'2. *wstak (istak) /okt 
'Vstak (istak) /ocm 


110 


2010 

3010 

3020 

3030 

3040 


check  the  number  of  initial  donors 
if (nidnr. eq.O) 

+  stop  '  BWACO:  ERROR  -  No  initial  donors  were  specified!' 
if (nidnr. It .nstak)  write (6, 2010) 
write (8,304fl)  nidnr 

return 

format  statements 

format (//'  BWACO:  WARNING  -  Some  stacks  are  not  initial  donors!') 
stak  called') 
region:  Ip5el0.3) 

stack;  ',4fl0.2) 
donors:  ',13) 


format  (/ 
format 
format 
format 

end 
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Figure  A- 16.  Subroutine  stak  firom  the  execution  module. 
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c 

c 

c 


c 

c 

c 

c 

c 

c 
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subroutine  blok  (ll:»nax,nstat,argiv,  lblok,nblok) 

determines  the  integer  block  length  and  nu^er 

of  blocks  that  yields  the  block  aspect  ratio 

that  most  closely  matches  a  given  aspect  ratio 

independent  variables: 

Ibmax  =  maximum  block  length 
nstat  =  number  of  stations 
argiv  »  given  aspect  ratio 

dependent  variables: 
iDlok  *  length  of  blocks 
nblok  =  number  of  blocks 

determine  the  minimum  block  length 
Ibmin^nst at /Ibmax 

initialize  the  counter 
numbr-0 

cycle  through  all  possible  block  lengths 
do  110  Iblck-lbmin, ibmax 

if  there  are  an  Integral  number  of  blocks  of  length  Iblck 
if (mod (nstat, Iblck) .eq.O)  then 

determine  the  number  of  blocks 
nblck=nst at /Iblck 


determine  the  block  aspect  ratio 
arblk=float  (Iblck) /float (nblck) 

as  soon  as  the  block  aspect  ratio 
exceeds  the  given  aspect  ratio 
if (arblk.ge. argiv)  then 

if  a  previous  solution  has  been  determined 
if (numBr.gt .0)  then 

compute  the  differences  between  the  given  aspect 
ratio  and  the  previous  and  current  aspect  ratios 
dold=abs (arglv-arblo) 
dcur=abs (argiv-arblk) 

if  the  given  aspect  ratio  is 
closer  to  the  previous  aspect  ratio 
if (dold.lt. dcur;  then 

corresponding 
pect  ratio 


return  the  block  length 
to  the  previous  block  a 
lblok*lblko 


otherwise 

else 

return  the  block  length  corresponding 
to  the  curent  block  aspect  ratio 
lblok=lblck 
nblok-nblck 


endlf 

ot herwlse 
else 

return  the  block  length  corresponding 
to  the  curent  block  aspect  ratio 
lblok=lblck 
nblok-nblck 


endlf 


return 


endlf 


Increment  the  counter  and 
numbr=numbr+l 
lblko=lblck 
nblko=nblck 
arblo-arblk 


save  the  current  values 


endlf 


continue 

If  a  previous  solution  exists 
If (numbr.gt .0)  then 

return  the  block  length  corresponding 
to  the  previous  block  aspect  ratio 
lblok=lblko 
nbl ok=nblko 
return 

otherwise 

else 

an  error  condition  exists  ,  j,, 

stop  '  BWACO:  ERROR  -  Region  cannot  be  blocked!' 


endlf 


end 
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Figure  A- 17.  Subroutine  blok  from  the  execution  module. 
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onooooooooooooooooooooooooooooo  oo 


subroutine  stns 
defines  the  stations 


station  variables: 
nstat  =  number  or  stations 
Iblok  =  length  of  station  bloks 
nblok  **  number  of  station  bloks 
xstat  *  array  of  x-coordinates  for  the  stations 

ystat  =  array  of  y-coordinates  for  the  stations 

zstat  =  z-coordinate  of  stations 

icoal  -  array  of  coalescence  statuses  for  the  stations 

wstat  *=  array  of  effective  yields  for  the  stations 

pstat  =  array  of  effective  peak  overpressures  for  the  stations 

region  variables: 
rinhb 


inhabited  build 
«  minimum  x-coord 
-  maximum  x-coord 

*  minimum  y-coord 

•  maximum  y-coord 


ng  radius 

nate 

nate 

nate 

nate 


xmln: 
xmax: 
ymln: 
ymax: 

local  variables: 
istan  -  station  index 
jstan  =  station  index 
Istat  “  station  counter 
xregn  *  x-dimension  of  region 
yregn  -  y-dimension  of  region 
asrat  «  aspect  ratio  of  region 
deltx  “  x-lncrement 
deity  »  y-increment 

common/conv/ 
common/icfg/  __ 
common/mxmn/  xm 


/conv/  okt , ocm, odsc. ocmkt3 
/icfg/  icgmb, iscrn# iprnt 
common/mxmn/  xmin.xmax, ymlh,ymax 
common/regn/  xcent, ycent, rlnnb,xmini, xmaxl,yroini, ymaxi 
common/ambi/  pambl,rambi 
common/sttn/  nstat 
common/sttb/  Iblok, nblok 
common/sttp/  xstat (1201, ystat (1201t2Stat 
common/sttv/  icoal (3600 J , pstat (3600) , wstat (3600) 


call 

e(8,3010) 


check  the  region  limits 
if (xmin.eg.xmax)  then 
xmin*xmini 
xmax«xmaxi 
endif 

if (ymln. eg. ymax)  then 

fin»yminl 
|x=ymaxt 

e (8, 3020)  xmin, xmax, ymin, ymax 

convert  the  region  limits  to  1-kt  standard  units 
xmln»xmin/ocm 
xmax*xmax/ocm 
ymin»»ymin/ocm 
ymax*ymax/ocm 

writel8,3020)  xmin, xmax, ymln, ymax 

define  the  number  of  stations 
nstat«36Q0 
lkmiax-120 

compute  the  coordinate  ranges  for  the  region 
xregn»xmax-xmin 
yregn-ymax-ymln 

distribute  the  station  Indices 
asrat =xre 
ca^  blok 
write  " 


olo^5l^ax?nstat, asrat ,  Iblok,  nblok) 
(8,3030)  l£>lok,r^^ -‘■ 


a  k./  aa  j.  a 

, nblok 

define  the  x  values 
deltx=xregn/f loat (Iblok) 
do  110  lstan«l, Iblok 

xstat  (Istan)  =  (float (istan-1) +0. 5) *deltx+xmln 
110  continue 

define  the  v  values 


120 


130 


3010 

3020 

3030 


define  the  y  values 
delty=yregn/tloat (nblok) 
do  120  3stan=l, nblok 

ystatt jstan)-  (float  ( jstan-1) +0. 5) *delty+ymin 
continue 

zstat=0. 

initialize  the  coalescence  status,  effective  yield, 
overpressure  and  ambient  pressure  and  density 
do  130  istat-l, nstat 
icoal  (lstat)=0 
wstat  (istat} -0. 
pstat (istat) -0. 
cohtinue 

call  at62(0.,  pambi, dummy, rambi, dummy) 
return 

format  statement 
format (/'  stns  called') 
format!'  region:  '.4fl0.2) 

format ('  blocks:'  ,5i4) 

end 
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Figure  A-18.  Subroutine  stns  from  the  execution  module. 
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c  _ _ _ _ _ 

c 

subroutine  rtim (ihbgn, imbgn, isbgn, ilbgn, ihend, imend, isend, ilend) 
c 

c  determines  and  reports  computation  time 

c 

c  independent  variables : 

c  ihbgn  =  hour  at  which  computation  begins 

c  imbgn  =  minute  at  which  computation  begins 

c  isbgn  =  second  at  which  computation  begins 

c  ilbgn  =  tenth  second  at  which  computation  begins 

c  ihend  =  hour  at  which  computation  ends 

c  imend  =  minute  at  which  computation  ends 

c  isend  =  second  at  which  computation  ends 

c  ilend  =  tenth  second  at  which  computation  ends 

c 

c  local  variables: 

c  tsbgn  =  time  at  which  computation  begins 

c  tsend  =  time  at  which  computation  ends 

c  tscmp  =  computation  time 

c 

c  compute  the  computation  time 

if ( ihend . It . ihbgn )  ihend=ihend+2 4 
tsbgn=3600 . *float (ihbgn) +60 . *float (imbgn) 

+  + (float (isbgn) +float (ilbgn) /lOO . ) 

tsend=3600 . *float (ihend) +60 .*float (imend) 

+  + (float (isend) +float (ilend) /lOO . ) 

t  scmp=t  send-t  sbgn 
c 

c  report  the  computation  time 

if (tscmp . It . 3600 . )  then 
if (tscmp . It . 60 . )  then 
write (6,2010)  tscmp 
else 

write ( 6, 2020)  tscmp/60. 
endif 
else 

write (6, 2030)  tscmp/3600. 
endif 
c 

return 

c 

c  format  statements 

2010  format (/'  Computation  Time  =',f7.2,'  seconds'/) 

2020  format (/'  Computation  Time  =',f7.2,'  minutes'/) 

2030  format  (/'  Computation  Time  =',f7.2,'  hours'/) 
c 

end 

c  _ _ _ _ _ 
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Figure  A- 19.  Subroutine  rtim  from  the  execution  module. 
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Subroutine  copt  is  shown  in  Figure  A-20.  It  cycles  through  all  possible  donors  (as  specified  by  user 
input)  computing  the  initiation  times  with  calls  to  prop  and  the  effective  yields  with  calls  to  effy.  It  also 
obtains  the  beginning  and  ending  times  of  the  computations  and  calls  rtim  to  compute  and  report  the  run 
time. 

PLOTTING  MODULE 

These  routines  plot  the  results  of  the  computations  on  the  request  of  the  user  using  proprietary 
graphics  routines. 

Subroutine  rest,  shown  in  Figure  A-21,  restores  the  original  units  and  eliminates  the  ground  plane 
factor. 

Subroutine  spmn,  shown  in  Figure  A-22,  displays  the  BWACO  plotting  menu  and  obtains  the  user’s 
request. 

Subroutine  pits,  which  is  not  listed,  generates  the  coalescence  and  peak  overpressure  maps  using  calls 
to  proprietary  routines. 

* 

interfac:e  module 

These  routines  set  up  the  session  and  direct  the  main  BWACO  functions  at  the  discretion  of  the  user. 

Subroutine  smmn,  shown  in  Figure  A-23,  displays  the  BWACO  main  menu  and  obtains  the  user’s 
request. 

Subroutine  exec,  shown  in  Figure  A-27,  calls  the  routines  of  the  execution  module  when  the  user 
requests  that  a  problem  be  run. 

Subroutine  rcfg,  shown  in  Figure  A-24,  reconfigures  the  session,  allowing  the  user  to  select  the  default 
display  and  printer. 
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subroutine  copt 

determines  the  worst-case  coalescence  pattern 

stack  variables: 
nstak  =  total  number  of  stacks 
idonr  =  array  of  initial  donor  switches 

local  variables: 
istak  =  stack  counter 

ihbgn  =  hour  at  which  computation  begins 
imbgn  =  minute  at  which  computation  begins 
isbgn  =  second  at  which  computation  begins 
ilbgn  =  tenth  second  at  which  computation  begins 
ihend  =  hour  at  which  computation  ends 
imend  =  minute  at  which  computation  ends 
isend  =  second  at  which  computation  ends 
ilend  =  tenth  second  at  which  computation  ends 

common/stkn/  nstak 
common/stkd/  idonr (100) ,nidnr 

log  call 
write  (8,3010) 

obtain  starting  time 
call  gettim( ihbgn, imbgn, isbgn, ilbgn) 

cycle  through  the  stacks 
idcnt=0 

do  110  istak=l, nstak 

if  the  stack  is  a  possible  initial  donor 
if (idonr (istak) .eq. 1)  then 

report  the  current  initial  donor  number 
idcnt=idcnt+l 
write  (6,2010)  idcnt,nidnr 

determine  the  propagation  sequence 
call  prop (istak) 

determine  the  (worst-case)  effective  yields 
call  effy 

endif 

110  continue 

obtain  completion  time 
call  gettim( ihend, imend, isend, ilend) 

determine  and  report  computation  time 

call  rtim (ihbgn, imbgn, isbgn, ilbgn,  ihend, imend, isend, ilend) 
return 


format  statements 
2010  format (  /'  Initial  Donor 
3010  format (/'  copt  called') 
c 

end 

c  _ 


',i3,'  of 


,i3,' :') 
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Figure  A-20.  Subroutine  coot  from  the  execution  module. 


74 


oo  o  ooo  oo  oo  oo  oooooooooooooooooooo  oo 


subroutine  rest 

restores  original  units  and  plots  the  results 

stack  variables : 
nstak  =  number  of  stacks 

xstak  =  array  of  x-coordinates  for  the  stacks 
ystak  =  array  of  y-coordinates  for  the  stacks 

station  variables: 

Iblok  =  length  of  station  blocks 
nblok  =  number  of  station  blocks 
nstat  =  number  of  stations 

wstat  =  array  of  effective  yields  for  the  stations 

local  variables: 
istak  =  stack  counter 
istan  =  station  index 
jstan  =  station  index 
istat  =  station  counter 


common/conv/ 
common/ stkn/ 
common/ stkp/ 
common/sttb/ 
common/sttn/ 
common/sttp/ 
common/ sttv/ 


okt, ocm, odsc, ocmkt3 
nstak 

xstak (100) , ystak (100) , zstak 

Iblok, nblok 

nstat 

xstat (120) , ystat (120) , zstat 
icoal (3600) ,pstat (3600) , wstat (3600) 


log  call 
write (8, 3010) 


restore  the  original  units  to  the  stack  coordinates 
do  110  istak=l, nstak 

xstak (istak) =ocm*xstak (istak) 
ystak (istak) =ocm*ystak (istak) 

110  continue 


restore  the  original  units  to  the  station  coordinates 
do  120  istan=l, Iblok 

xstat (istan) =ocm*xstat (istan) 

120  continue 

do  130  jstan=l, nblok 

ystat ( jstan) =ocm*ystat (jstan) 

130  continue 


restore  the  original  units  to  and  eliminate  the 
ground  plane  factor  from  the  effective  yields 
do  140  istat=l, nstat 

wstat (istat) =0 . 5*okt*wstat (istat) 

140  continue 


return 


format  statement 
3010  format (/'  rest  called') 
c 

end 


c 
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Figure  A-21.  Subroutine  rest  from  the  plotting  module. 
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110 
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subroutine  spmn (irqst, iptyp, idevi) 

displays  plotting  menu  and  obtains  user  request 

dependent  variables : 
irqst  =  user  request  index 
iptyp  =  index  of  requested  plot  type 
idevi  =  index  of  requested  device 

log  call 
write (8, 3010) 

display  menu 
write (6,2010) 
write (6,2020) 
write (6, 2030) 
write (6,2040) 
write (6,2050) 
write (6,2060) 
write (6,2070) 
write (6,2080) 

obtain  user  request 
read(5, 1010,err=110)  irqst 
if (irqst. It. 0. or. irqst. gt. 5)  go  to  110 
write (8, 3020)  irqst 

assign  requested  plot  type 
if (irqst .eq. 1 .or . irqst .eq.2)  then 
iptyp=2 
else 

iptyp=l 

endlf 

assign  requested  device 
if (irqst . eq . 1 . or . irqst . eq . 3 )  then 
idevi=iscrn 
else 

idevi=iprnt 

endif 


return 


format  statements 
format (il) 
format (// 


format ( 
format  ( 
format ( 
format ( 
format  ( 
format  ( 
format ( 
format (/ 
format (' 


/' 


BWACO 

Enter 


'\) 

spmn  called' 
irqst=' il) 


PLOTTING  MENU') 

1  to  plot  Coalescence  Map  (Screen)') 

2  to  plot  Coalescence  Map  (Printer) ' ) 

3  to  plot  Pea]c  Overpressure  Map  (Screen) ' ) 

4  to  plot  Pea)c  Overpressure  Map  (Printer) ' ) 

5  to  reconfigure  the  session' ) 

0  to  return  to  the  Main  Menu' ) 


) 


end 
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Figure  A-22.  Subroutine  spmn  from  the  plotting  module. 
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subroutine  sininn  (iscur,  irqst) 

displays  main  menu  and  obtains  user  request 

display  menu  heading 
write (6,2010) 


check  for  presence  of  bwaco.dat 
open (7, f ile='bwaco .dat' , status=' old' , err=120) 
close  (7) 

display  complete  menu 
write (6,2020) 
if (iscur .eq. 0)  then 
write (6,2030) 

0lS0 

write (6,2031) 
endif 

write (6,2041) 
write (6,2050) 
write (6,2060) 

obtain  user  request 
read(5, 1010,err=110)  irqst 

check  for  valid  request 
if (irqst . It . 0 .or . irqst .gt . 3)  go  to  110 

go  to  130 

display  partial  menu 
write (6,2020) 
write  (6, 2042) 
write (6, 2050) 
write  (6, 2060) 

obtain  user  request 
read(5, 1010,err=120)  irqst 

check  for  valid  request 
if (irqst . It . 0 .or . irqst .gt .2)  go  to  120 

adjust  request  number 
if (irqst .eq. 2)  irqst=3 

return 


format  statements 
format (il) 

format  (//'  BWACO  MAIN  MENU  ') 
format  (/'  Enter  1  to  run  the  current 
format  (  '  2  to 

format  (  '  2  to 

format  (  '  3  to 

format  (  '  2  to 

format  (  '  0  to 

format  (  '  '\) 


problem' ) 
plot  the  previous  problem' ) 
plot  the  current  problem' ) 
reconfigure  the  session' ) 
reconfigure  the  session' ) 
end  the  session' ) 


end 
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Figure  A-23.  Subroutine  smmn  Irom  the  interface  module. 
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if(letuif.n«.0)  fo  110 


ills 
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7J  format  #  to  Ktaro  to  tM  Mno  » 

ie73  format  4  to  rotsm  to  tbo  flottlay  ■odd') 

010  format  ' 


isss  !ss:s 
}i;i  ss;::s 
ISii  iiXi 
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tlfO^fopm^ 

S13e^|o|fm^ 
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im  fes:u'/«‘fct:?t8:*s.:’4?s»jr/ii^:  -m 
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**•*♦5*^^:  W  Think  Jot  Lo  Dona  t>:  IT  think  Jat  UDmo') 

**'’*4|*F*|i;  n  Falat  Jot  io  Dona  t7:  W  Point  Jot  U  Dona') 

Slie^fof^,  Portaorlpt  Printor  ft:  toofoonlotoJ  Pootooriyt') 

m  KSRjf'  R>r«S!orT  ~ 


1«40  format  j'  noofiftrotloii  • 


OO  OO  O  O  OO  O  OOO  OOO  O  OO  OO  O  O  OO  O  OOO  00 


subroutine  savn 
saves  new  data  to  file 


common/icfg/ 
common/totl/ 
common/ulab/ 
common/stkn/ 
common/ stkp/ 
common/stkd/ 
common/regn/ 
common/sttn/ 
common/sttb/ 
common/ sttp/ 
common/sttv/ 
common/mxmn/ 

character*80  ulabl 

log  call 
write (8, 3010) 


icomb, iscrn, iprnt 

wtotl 

ulabl 

nstak 

xstak  (100)  ,ystak  (100) ,  zstak 
idonr (100) , nidnr  .  .  ... 

xcent,ycent, rinhb, xmini, xmaxi,ymini, ymaxi 
nstat 

Iblok.nblok  _ _ 

xstat  (120l,ystat  (120)  ,  zstat 
icoal  (slooJlpstat  (36o6)  ,wstat  (3600) 
xmi n , xmax , ymi n ,  ymax 


open {7. file='bwaco.dat' , form=' unformatted' ,err=110) 
write  (7)  ulabl . 


write 

write 

write 

write 

write 

write 

write 

write 

write 

write 

write 

close 

return 


nstat, Iblok, nblok, nstak, nidnr 
wt  o 1 1 , xmi n , xmax,  ymi n , ymax 
xcent, ycent, rinhb 
'xstak  (i) ,  i=l,  nstak 
ystak ' i) , i=l, nstak 
idonr ii),i=l, nstak 
xstat I i) , i=l, Iblok 
ystat I i ) , i=l, nblok 
pstat I i ) , i=l, nstat 
wstat I i) , i=l, nstat 
icoal (i) , i=l, nstat 


error  detected  ,  ,  ,  ,  . , 

110  stop  '  BWACO:  ERROR  -  Unable  to  open  bwaco.dat!' 

format  statement 
3010  format (/'  savn  called') 

end 


subroutine  geto 

gets  old  plotting  data  from  file 

common/icfg/  icomb, iscrn, iprnt 
common/totl/  wtotl 
common/ulab/  ulabl 

common/stkn/  nstak  _ _  ,  _ 

common/stkp/  xstak  100) , ystak (100) , zstak 
common/stkd/  idonr (100) , nidnr  .  .  ...  , 

common/regn/  xcent, ycent, rinhb, xmini, xmaxi, ymini, ymaxi 
common/sttn/  nstat 
common/sttb/  Iblok, nblok 
common/sttp/  xstat  120), ystat (120), zstat 
common/sttv/  icoal (3600) , pstat (3600) , wstat (3600) 
common/mxmn/  xmin, xmax,  y min, ymax 

character*80  ulabl 

log  call 
write (8, 3010) 


open 

read 

read 

read 

read 

read 

read 

read 

read 

read 

read 

read 

read 

close 


7 , end= 
7, end= 
7,  end- 
7,  end 
7,  end 
7,  end 
7,  end 
7,  end 
7,  end 
7,  end 
7,  end 
7.  end 
{^) 


bwaco 
130, err 
130, err 
130, err 
130, err 
130  err 
130;err 
130, err 
130, err 
130, err 
130, err 
130;err 
130, err 


120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 

120 


ulabl 

nstat, Iblok, nblok, nstak,  nidnr 
wtotl, xmin , xmax,  ymi n , y max 
xcent, ycenp , rinhb 


xstak 
I  ystak 
idonr 
xstat 
ystat 
pstat 
wstat 
icoal 


, 1=1, nstak 
, i=l, nstak 

5| 

, 1=1, nblok 
, i=l, nstat 
, i=l, nstat 
, i=l, nstat 


return 


errors  detected  ,  ,  ,  ,  , . 

110  stop  '  BWACO:  ERROR  -  Unable  to  open  bwaco.dat  ' 

120  stop  '  BWACO:  ERROR  -  Unable  to  read  bwaco.dat!'.. 

130  stop  '  BWACO:  ERROR  -  Premature  end  of  bwaco.dat!' 

format  statement,.  ,  , 

3010  format!/'  geto  called') 

end 
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Figure  A-25.  Subroutines  savn  and  geto  from  the  interface  module. 
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Subroutine  savn,  shown  in  Figure  A-25,  saves  the  results  of  a  computation  to  a  file. 

Subroutine  geto,  shown  in  Figure  A-25,  gets  the  results  of  a  previously  run  computation  from  a  file. 

Subroutine  pmen,  shown  in  Figure  A-26,  processes  requests  from  the  plotting  menu  when  a 
computation  is  complete  or  when  the  user  requests  plotting  from  the  main  menu. 

Subroutine  cnfg,  shown  in  Figure  A-27,  obtains  the  default  configuration  parameters  from  the  file 
bwacoxfg  or  calls  rcfg  if  bwaco.cfg  does  not  exist. 

Subroutine  cnst,  shown  in  Figure  A-28,  sets  the  required  constants. 

Subroutine  mmen,  shown  in  Figure  A-26,  processes  requests  fiom  the  main  menu. 

Program  bwac,  shown  in  Figure  A-27,  is  the  main  program.  It  configures  the  session  with  calls  to 
crtfg  and  cnst,  and  initiates  the  main  menu. 
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subroutine  pmen 

solicits  and  processes  a  user  request  from  the  plotting  menu 

local  variables: 
irqst  =  user  request  index 
iptyp  =  index  of  requested  plot  type 
idevi  =  index  of  requested  device 

log  call 
write  (8,3010) 

restore  original  units 
call  rest 


110 


sj^icit  user  request 


spmn (irqst, iptyp,  idevi) 

if (irqst . gt . 0 . and. irqst . It .5) 

process  plotting  request 
call  pits (iptyp, idevi ) 

go  to  110 

endif 


then 


if (irqst .eq. 5)  then 

process  reconfiguration  request 
call  rcfg(0,3) 
call  savn 


go  to  110 
endif 


return 


3010 


format  statement 
format (/'  pmen  called') 


end 


subroutine  mmen 


solicits  and  processes  a  user  request  from  the  main  menu 

local  variable: 
irqst  =  user  request  index 

log  call 
wriee(8,3010) 

solicit  user  request 
iscur=0 

110  call  smmn (iscur, irqst) 

if (irqst .eq. 1)  then 

mcess  a  request  to  run  a  problem 
exec 
call  savn 
call  pmen 
iscur=l 

endif 


if (irqst .eq. 2)  then 


process  a  request  to  plot  a  problem 
call  geto 
gall  pmen 
iscur=l 


endif 


if (irqst. eq. 3)  then 

process  a  request  for  session  reconfiguration 
call  rcfg(0,2T 

endif 


if (irqst . ne . 0)  go  to  110 
return 


c 

c 


format  statement 
format (/'  mmen  called') 

end 


P000090 

POOOlOO 

POOOllO 

PQ0012Q 

P000130 

P000140 


mm 


0230 


160 


P000320 

P000330 


P000380 

P000390 

PQ0040g 

PQg0410 

P000420 

P000430 


M06062g 

Mgooo3o 

Mggoo40 


M000080 


Mgo 

MOO 

mBo 

MOO 

MOO 


30186 

o6190 

00200 

30210 

10220 

im 

)250 

260 

5270 


MooSiiS 

M000340 


0390 

34gg 


M000450 

Mg00460 

M000470 

Mgoo48g 

MQogpo 

mOoOsoo 

M000510 


Figure  A-26.  Subroutines  pmen  and  mmen  from  the  interface  module. 
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110 


1010 

3010 


program  bwac 

is  the  BWACO  main  program 

report  the  BWACO  version 
wrlte(6,2010) 

open  the  session  log  _ _ 

open (8, file-' bwaco.log' , form-' formatted' , err=110) 
write (1,3010) 

configure  the  session 
call  cnfg 

define  the  constants 
call  cnst 

run  the  main  menu 
call  mmen 

close  the  session  log 
write (8,3020) 
close (8} 

terminate  execution 
stop  '  BWACO:  Session  Terminated.' 

error  detected 

stop  '  BWACO:  ERROR  -  Unable  to  open  bwaco.log!' 

format  statement 
format (//'  BWACO:  May  1994') 
formate'  BWACO  Session  Log') 
format!/'  termination  requested') 

end 


110 


2010 

3010 

3020 


subroutine  cnfg 

configures  the  session 

configuration  variables; 
icomb  “  pressure  combination  algorithm  switch 
iscrn  *  default  display  switch 
iprnt  =  default  printer  switch 

common/lcfg/  icomb, iscrn, iprnt 

log  call 
write (8, 3010) 

read  data  from  the  configuration  file 
open  (7,  flle»'*bwaco,cfg*',  form^^formatted**,  status»**old",err-110) 
read (7l 1010,end*110,err=110)  icomb, iscrn, iprnt 
close  (7) 
go  to  120 

generate  default  configuration  file 
icomb*2 
iscrn*10 
iprnt“28 
call  rcfg(l,l) 


120  return 


format  statements 
format (11,212) 
format!/'  cnfg  called') 

end 


subroutine  exec 

sets  up  and  executes  the  coalescence  computation 


log  call 
write (8,3010) 

read  problem  input 
call  inpt 

define  the  stacks 
call  stak 

define  the  stations 
call  stns 

make  the  coalescence  computation 
call  copt 

return 


format  statement 
3010  format!/'  exec  called') 

end 
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Figure  A-27.  Program  hwac  and  subroutines  cnf2  and  exec  from  the  interface  module. 
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3010 

c 

c 


subroutine  cnst 

defihes  the  constants 

common/conf /  thrd, tpi 
common/conv/  okt, ocm, odsc, ocmktS 
common/coni/  gamma , gmone , gpone 
common/cona/  pl,cl,rl,tl 
common/conp/  a0p,b0p2, cOp, alp,blp, rOp 
common/cons/  r0s,c0s,d0s 
common/conk/  ak0,ak(o) 

log  call 
write (8, 3010) 

fixed  constants 
thrd=l . /3 . 
tpi=8 . *atan (1 . ) 
twthrd=2 . **thrd 

conversion  factor  constants 
okt=2 .e6 
ocm=0 . 01 
odsc=68947 . 6 

ocmkt3=2 . 604004e-4*twthrd 

ideal  gas  constants 
gamma=l . 404574 
gmone=gamma-l . 
gpone=gamma+l . 

sea-level  atmospheric  constants 
call  at62(0.,  pl,cl,rl,tl) 

positive  phase  duration  constants 

(time  in  s,  distance  in  cm,  weight  in  kt) 
aOp  =4,313976E-05 
b0p2=6.261821E+08 
cOp  =3,350000E-01 
alp  =5.361303E+12 
blp  =3.507575E+00 
rOp  =4,490000E+03 

shock  arrival  time  constants 

(time  in  s,  distance  in  cm,  weight  in  kt) 
r0s=l . 6195e4 
c0s=3.3897e4 
d0s=1.5000e4 

kingery  overpressure  fit  constants 

(pressure  in  lb/in**2,  time  in  s,  distance  in  ft,  weight  in  lb) 
akO  =  7.0452041 
ak(l)=-l. 6277561 
ak(2)=-0. 27399088 
ak(3)=-6.5973136e-2 
ak(4)=  6.5412563e-3 
ak(5)=  4.8236359e-2 
ak(6)=-2.0072553e-2 
ak(7)=  3.0190449e-3 
ak(8)=-1.5984026e-4 

return 

format  statement 
format (/'  cnst  called') 

end 
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Figure  A-28.  Subroutine  cnst  from  the  interface  module. 
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INTRODUCTION 


This  manual  provides  background  information  on  the  BWACO  methodology, 
outlines  hardware  and  software  requirements  and  directs  installation  of  the  code.  The 
BWACO  Tutorial  gives  instructions  for  running  sample  problems  and  creating  and  running 
your  own  problems.  In  addition,  the  BWACO  Reference  describes  all  BWACO  functions 
and  features  in  alphabetical  order.  Finally,  a  detailed  description  of  BWACO,  listing  most 
of  the  source  code,  is  provided.  Users  who  experience  difficulty  should  contact  the 
authors  at  (410)  278-6214/6241. 


BACKGROUND 

The  BWACO  program  was  written  to  predict  regions  of  blast  wave  coalescence 
associated  with  the  sequential  mass  detonation  of  distributed  ammunition  stacks 
positioned  as  specified  by  the  user. 

It  accomplishes  this  by  applying  the  assumption  that  communication  of  mass 
detonation  between  stacks  is  controlled  by  a  user-supplied  propagation  velocity  (which 
might  represent  the  velocity  of  fragments  from  the  donor  stack).  BWACO  uses  this 
propagation  velocity  to  determine  the  order  and  timing  of  detonation  of  the  stacks.  It 
recognizes  that  the  identity  of  the  initial  donor  is  not  generally  known  and,  unless  a  single 
initial  donor  has  been  specified  by  the  user,  considers  all  possible  initial  donors  in  order 
to  predict  the  worst-case  environment  at  any  point. 

Then,  under  the  assumption  that  blast  waves  propagate  independently  of  one 
another,  BWACO  applies  the  criterion  that  coalescence  occurs  wherever  and  whenever 
the  shock  associated  with  one  blast  wave  encroaches  into  the  positive  overpressure 
phase  of  another  blast  wave.  This  assumption  generally  causes  premature  identification 
of  coalescence.  When  coalescence  is  detected  at  a  point,  it  determines  the  combined 
peak  overpressure  by  first  combining,  at  their  center  of  charge,  the  stacks  producing  the 
waves  which  coalesce  and  then  computing  the  peak  overpressure  associated  with  the 
combined  stacks  at  that  point.  This  method  is  suitable,  at  least,  for  closely  spaced 
stacks. 


BWACO  graphically  reports  regions  of  coalescence,  showing  the  effective  total 
explosive  weight  and  the  peak  overpressure  throughout  the  region  of  interest. 

In  applying  BWACO,  it  is  important  to  remain  aware  of  limitations  imposed  by  the 
assumptions.  The  accuracy  of  BWACO  has  not  been  assessed  except  in  the  case  of 
closely  spaced  charges.  As  the  separation  between  stacks  is  increased,  the  accuracy, 
especially  that  of  the  peak  overpressure  predictions,  is  expected  to  decrease. 

Detailed  background  information  can  be  found  in  the  references. 
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HARDWARE  AND  SOFTWARE  REQUIREMENTS 


BWACO  should  run  in  the  DOS  environment  on  any  personal  computer  with  an 
80786  microprocessor.  However,  to  minimize  the  running  times  for  large  problems,  a 
high-speed  computer  with  at  least  an  80486  microprocessor  is  recommended. 

An  editor  (such  as  the  DOS  editor)  that  produces  pure  ASCII  text  (without  hidden 
characters)  is  required  to  create  problem  input  files.  A  word  processing  program  shouid 
not  be  used  unless  the  file  can  be  saved  in  ASCII  format. 

BWACO’s  graphics  supports  a  number  of  different  video  displays  and  printers.  The 
displays  are  listed  in  the  following  table: 


Type 

Mode 

Resolution 

Colors 

CGA 

4 

320x200 

4 

CGA 

5 

320x200 

4 

CGA 

6 

640x200 

BW 

Hercules 

720x348 

BW 

EGA 

13 

320x200 

16 

EGA 

14 

640x200 

16 

EGA 

15 

640x350 

BW 

EGA 

16 

640x350 

16 

VGA 

17 

640x480 

BW 

VGA 

18 

640x480 

16 

ATI 

81 

640x480 

16 

ATI 

82 

752x410 

16 

ATI 

83 

800x560 

16 

The  following  printers  are  supported: 

single-density  printer 
double-density  printer 
quad-density  printer 
Epson  printer  (plot  mode) 

Postscript  printer 
Hewlett-Packard  Laser  Jet 
Hewlett-Packard  Think  Jet 

Hard  copies  requested  during  a  BWACO  session  are  automatically  sent  to  the  printer  at 
the  end  of  the  session. 
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INSTALLATION 


Please  follow  these  installation  instructions  very  carefully.  Before  installing 
BWACO,  create  a  directory  to  hold  the  BWACO  files.  For  example,  you  can  create  the 
directory  bwaco  by  entering  the  following  line: 

>mkdir  bwaco 

Make  that  directory  the  current  directory. 

>cd  bwaco 

Insert  the  BWACO  installation  diskette  into  the  appropriate  drive.  If  you  inserted 
it  into  drive  a:,  enter  the  following  line: 

>a: install  a 

If  you  inserted  it  into  drive  b:,  enter  the  following  line: 

>b: install  b 

This  will  cause  four  files  to  be  copied  to  your  new  directory,  after  which  the  installation 
is  complete. 


BWACO  TUTORIAL 


Getting  Started 

BWACO  always  looks  for  input  in  a  file  with  the  extension  inp.  A  sample  input 
file  has  been  supplied  and  running  this  problem  before  trying  your  own  is  recommended. 
The  file  three. inp  describes  three  25,000-lb  stacks  of  ammunition,  separated  by  200 
meters,  lying  at  the  vertices  of  an  equilateral  triangle.  To  start  a  BWACO  session  and 
run  this  problem,  simply  enter  the  following  line: 

>bc  three 

The  first  time  you  run  BWACO  the  Configuration  Menu  will  appear.  This  menu, 
which  is  similar  to  the  other  BWACO  menus,  requires  you  to  enter  a  single-digit  number 
to  select  one  of  the  menu  options: 
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BWACO  CONFIGURATION  MENU 

Default  Display:  VGA  Mode  18  640x480/16 

Default  Printer:  Postscript  Printer 

Enter  1  to  change  the  default  display 
2  to  change  the  default  printer 
0  to  proceed  to  the  Main  Menu 

The  Configuration  Menu  allows  you  to  tell  BWACO  which  display  and  printer  are  installed 
on  your  system.  The  default  display  is  vga  Mode  18  and  the  default  printer  is 
Postscript  Printer.  If  your  display  or  printer  is  different,  make  changes  now.  If  you 
do  not  specify  a  display  or  printer,  BWACO  will  save  and  use  these  defaults.  To  change 
the  display  or  printer,  simply  enter  either  l  or  2  and  make  your  selection  from  the  list 
displayed.  BWACO  will  remember  your  specifications  in  all  subsequent  sessions.  To 
leave  the  menu,  simply  strike  the  enter  key.  Typing  0  is  not  necessary.  After  you  have 
completed  the  initial  configuration,  the  Main  Menu  will  appear: 

BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 
2  to  reconfigure  the  session 
0  to  end  the  session 


To  run  the  problem  configuration  described  in  three .  inp,  enter  1.  You  will  have 
to  wait  some  period  of  time  while  the  problem  runs.  BWACO  will  keep  you  informed  of 
what  it  is  doing.  When  the  computation  is  complete,  BWACO  will  tell  you  how  long  it  took 
to  run  on  your  system  and  the  Plotting  Menu  will  appear: 

BWACO  PLOTTING  MENU 

Enter  1  to  plot  the  Coalescence  Map  (Screen) 

2  to  plot  the  Coalescence  Map  (Printer) 

3  to  plot  the  Peak  Overpressure  Map  (Screen) 

4  to  plot  the  Peak  Overpressure  Map  (Printer) 

5  to  reconfigure  the  session 

0  to  return  to  the  Main  Menu 


You  can  read  about  the  Coalescence  and  Peak  Overpressure  Maps  in  the 
reference  section  of  this  manual.  Produce  the  Coalescence  Map  on  your  screen  by 
entering  l.  Strike  any  key  to  return  to  the  Plotting  Menu.  Enter  2  to  obtain  a  hard  copy. 
All  hard  copies  requested  during  a  session  will  be  printed  when  the  session  is  terminated. 
Try  the  Peak  Overpressure  Map  in  a  similar  fashion.  If  you  experience  problems  with 
plotting,  you  may  have  made  the  wrong  configuration  selections.  You  can  correct  them 
by  entering  5  from  the  Plotting  Menu. 
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When  you  are  finished  plotting,  return  to  the  Main  Menu: 

BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 

2  to  plot  the  current  problem 

3  to  reconfigure  the  session 
0  to  end  the  session 

This  menu  now  provides  the  additional  option  of  plotting  the  results  of  the  problem  without 
rerunning.  End  the  session  by  striking  the  enter  key.  You  may  receive  a  query 
regarding  your  printer.  If  so,  simply  strike  enter  again.  Any  hard  copies  you  requested 
will  print  at  this  time. 


Creating  Your  Own  Input  Files 

BWACO  requires  its  input  file  to  have  a  very  specific  format.  Before  creating  a 
new  input  file,  let’s  look  at  the  format  of  the  file  three .  inp. 

Three  25-klb  Stacks,  200-m  Separation  (three. inp) 

1  — Xmin —  1  I  — Xmax —  |  |  — Ymin —  |  |  — Ymax —  I 

-700.  700.  -700.  700. 

I - Xs - II - Ys - II - Ws - II - ^Vs - 1  D 

100.  -57.74  25000.  2500.0  1 

-100.  -57.74  25000.  2500.0  1 

0.  115.48  25000.  2500.0  1 

The  first  line  of  this  file  contains  a  descriptive  problem  title  of  the  user’s  choice.  It  may 
contain  up  to  80  characters.  It  is  used  to  label  the  plots.  As  a  convenience,  we  have 
included  the  input  file  name  in  parentheses  at  the  end  of  the  title.  The  second  line  is  a 
ruler  line  delimiting  the  fields  in  which  the  bounds  of  the  region  of  interest  are  placed. 
These  bounds,  which  appear  on  the  third  line,  determine  the  region  in  which  the  stations 
are  distributed.  The  fourth  line  is  another  ruler  line  specifying  the  fields  in  which  stack 
data  is  placed.  The  three  stacks  are  described,  one  to  a  line,  on  the  remaining  lines. 
The  stack  coordinates,  the  TNT  equivalent  explosive  weight,  and  the  propagation  velocity 
are  specified.  In  addition,  the  D  (initial  donor)  switch  is  set  to  1  to  ensure  that  each  stack 
may  act  as  an  initial  donor. 

Let’s  see  how  problem  input  is  generated.  The  net  explosive  weight  for  each  stack 
should  be  expressed  in  equivalent  pounds  of  TNT,  distances  in  meters  and  propagation 
velocities  in  meters  per  second.  Suppose  we  have  four  25,000-lb  ammunition  stacks 
arranged  as  shown  in  Figure  B-1.  First,  we  must  define  a  Cartesian  coordinate  system. 
The  orientation  and  origin  location  of  our  coordinate  system  is  arbitrary  and  may  be 
chosen  to  suit  our  convenience.  One  possibility,  along  with  the  coordinates  of  each  stack, 
is  shown  in  Figure  B-2.  In  addition  to  this,  we  also  have  to  know  the  propagation  velocity 
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associated  with  each  stack.  Since  no  better  guidance  is  presently  available,  we  have 
fixed  this  at  2,500  m/s. 

We  now  know  enough  to  generate  the  input  file.  In  order  to  create  your  own 
problem  input  files,  you  will  need  to  use  an  editor  that  produces  pure  ASCII  files  (without 
hidden  characters),  such  as  the  DOS  editor.  A  word  processing  program  should  not  be 
used.  As  a  convenience,  an  input  template  has  been  included  with  BWACO.  To  start 
out,  simply  copy  this  file  to  the  new  input  file,  as  follows: 

>copy  template . inp  four.inp 

Then,  edit  the  file  four.inp.  Initially,  it  looks  like  this: 

I  — Xmin —  I  I  — Xmax —  1  1  — Ymin —  |  1  — Ymax —  | 

I - Xs - 1  I - Ys - 1  1 - Ws - 1  I - Vs - 1  D 

First  we’ll  add  a  descriptive  title: 

Four  25-klb  Stacks,  200-m  Separation  (four.inp) 

I  — Xmin —  I  I  — Xmax —  1  |  — Ymin —  |  |  — Ymax —  | 

I - Xs - 1  1 - Ys - 1  I - Ws - 1  1 - ^Vs—  I  D 

Next,  since  we  don’t  know  how  big  the  region  of  interest  should  be,  we’ll  let  BWACO 
determine  its  limits  by  specifying  X  and  Y  bounds  which  are  equal: 

Four  25-klb  Stacks,  200-m  Separation  (four.inp) 

I  — Xmin —  I  I  — Xmax —  |  |  — Ymin —  |  |  — Ymax —  | 

0.  0.  0.  0. 

I - Xs - 1  I - Ys - 1  I - Ws - 1  1 - ^Vs - 1  D 

Finally,  we’ll  specify  the  stack  coordinates,  weights  and  propagation  velocities,  setting  the 
initial  donor  switch  to  1  for  each  stack: 


Four  25-klb  Stacks,  200-m  Separation  (four.inp) 
I  — Xmin —  I  I  — Xmax —  |  |  — Ymin —  |  I  — Ymax —  | 

0.  0.  0.  0. 


-100. 

100. 

-100. 

100. 


-100. 

-100. 

100. 

100. 


- Ws - 1  I - ^Vs - 1 

25000.  2500. 
25000.  2500. 
25000.  2500. 
25000.  2500. 


D 

1 

1 

1 

1 


This  completes  the  generation  of  the  input  file  four.  inp.  Try  running  this  problem: 


>bc  four 
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A 


A 


A  A 


Figure  B-1.  Arrangement  of  four  25,000-lb  ammunition  stacks. 


y 


(•100.100)^  ^(100.100) 
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(•lOO.'lOO)^  ^(100.-K)0) 


Figure  B-2,  A  coordinate  system  for  this  arrangement. 
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The  plots  produced  will  be  somewhat  distorted.  This  occurs  because  the  algorithm  used 
to  extend  the  axes  to  a  natural  termination  point  treats  the  X  and  Y  axes  differently.  To 
eliminate  the  distortion,  it  is  necessary  to  remn  the  problem  after  modifying  the  input  file 
to  specify  region  limits: 


Four  25-klb  Stacks,  200-m  Separation  (four.inp) 
I  — Xmin —  I  I  — Xmax —  |  |  — Ymin —  |  1  — Ymax —  | 

-800.  800.  -800.  800. 


I - Xs - 1  I - Ys - 1 

-100.  -100. 

100.  -100. 

-100.  100. 

100.  100. 


- Ws - 1  1 - ^Vs - 1 

25000.  2500. 
25000.  2500. 
25000.  2500. 
25000.  2500. 


D 

1 

1 

1 

1 


With  this  modification  (which  is  not  always  necessary  when  allowing  BWACO  to 
determine  the  extent  of  the  region  of  interest)  undistorted  plots,  as  shown  in  Figures  B-3 
and  B-4,  are  produced. 


Obtaining  Piots  From  Eariier  Computations 

We  can  easily  obtain  new  plots  from  the  problem  described  in  three. inp  run 
earlier.  Be  sure  to  wait  until  all  previously  requested  plots  have  printed,  otherwise  the 
files  may  be  deleted  before  they  print.  Simply  enter  the  following  line: 

>bc  three 

The  Main  Menu  appears  as  follows: 

BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 

2  to  plot  the  previous  problem 

3  to  reconfigure  the  session 
0  to  end  the  session 

The  previously  obtained  results  may  be  plotted  by  entering  2  and  making  the  desired 
selections  from  the  Plotting  Menu.  End  the  session  by  striking  enter  twice. 
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Figure  B-3.  Coalescence  Map. 
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100.0 

100.0 

25.0 

0.0 

0.0 

100.0 

Blast  Hazard  Prea  (km*) 


-  RR  3B5-64:  1 . 5055 
BURCO:  0.6050 
Contour  Levels  (psi): 


1.2  1.7  2.3 
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U.S.  Rryesrch  Laboratory 

BNflCO 

Peak  Overpressure  Map 

06^06/1994  16;00 


Figure  B-4.  Peak  Overpressure  Map. 


96 


BWACO  REFERENCE 


Coalescence  Map 

An  example  of  a  Coalescence  Map  is  shown  in  Figure  B-3.  Regions  in  which 
coalescence  is  detected  are  hatched  and  explosive  weight  contours  are  drawn.  Each 
contour  is  labeled  to  indicate  the  total  explosive  weight  contributing  to  the  worst-case 
blast  environment  within  that  contour.  The  outermost  contour  follows  the  0.9  or  1 .2  psi 
overpressure  limit  for  inhabited  buildings  defined  in  Army  Regulation  (AR)  385-64.  For 
comparison  purposes,  a  circle,  centered  at  the  center  of  charge  of  all  the  stacks  and 
having  the  inhabited  building  distance  for  the  combined  explosive  weight  of  all  the  stacks 
as  its  radius,  is  also  plotted.  As  many  contour  levels  as  space  permits  are  listed  in  the 
legend.  (See  Peak  Overpressure  Map  and  Plotting  Distortion.) 


Contour  Levels 

The  contour  levels  used  in  the  Coalescence  Map  are  listed  in  the  legend  (space 
permitting).  These  levels  describe  the  total  explosive  weight  contributing  to  the  worst- 
case  blast  environment  within  the  contour.  The  contour  levels  used  in  the  Peak 
Overpressure  Map  are  those  specified  in  AR  385-64  for  the  protection  of  a  variety  of 
assets. 


Files  and  Program  Execution 

The  file  be. bat  is  the  DOS  batch  file  that  executes  BWACO  and  manages  its 
files.  The  file  bwaco.exe  is  the  executable  BWACO  program.  The  files  * .  inp  are  the 
user’s  problem  specification  input  files  (where  *  represents  the  user’s  descriptive  file 
name).  (See  Input  Files.)  The  files  *  .dat  contain  the  plotting  data,  enabling  BWACO 
to  plot  results  without  rerunning  problems.  The  files  *.log  contain  logs  of  BWACO’s 
most  recent  session  pertaining  to  the  problem  identified  by  *.  The  user  starts  a  BWACO 
session  by  entering  be  *.  The  batch  file  first  deletes  any  old  plotting  files.  For  this 
reason,  a  BWACO  session  should  never  be  started  until  all  plots  from  the  previous 
session  have  printed.  Next,  it  copies  * . inp  and  * . dat  (if  it  exists)  to  bwaco. inp  and 
bwaco.dat,  respectively.  It  then  executes  bwaco.exe,  which  creates  the  files 
bwaco.dat  and  bwaco.log.  On  termination  of  the  session,  the  batch  file  deletes 
bwaco .  inp,  moves  bwaco .  dat  and  bwaco .  log  to  * .  dat  and  * .  log,  respectively 
and  sends  all  plotting  files  to  the  printer  using  the  print  command. 
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Inhabited  Building  Limit 

AR  385-64  defines  permissible  overpressure  exposure  levels  for  a  variety  of 
assets.  The  lowest  such  levels  are  those  commonly  associated  with  inhabited  buildings 
(0.9  to  1.2  psi).  The  levels  are  also  given  as  scaled  distances  using  classical  blast 
scaling  laws.  The  minimum  distance  is  further  limited  by  consideration  of  hazardous 
fragment  density.  BWACO  combines  all  the  stacks  by  summing  their  explosive  weights 
and  determining  their  center  of  charge.  Using  the  combined  weight,  it  computes  the 
inhabited  building  distance.  Using  the  center  of  charge,  it  plots  the  combined  inhabited 
building  circle  on  the  Coalescence  Map.  On  the  Peak  Overpressure  Map,  it  plots  a 
contour  produced  by  centering  inhabited  building  distance  circles  at  each  stack.  This 
represents  the  method  currently  required  by  the  regulation  and  allows  comparison  with 
BWACO’s  predictions. 


Initial  Donors 

When  detonation  propagates  from  one  stack  to  another,  the  first  stack  of  such  a 
pair  is  referred  to  as  the  donor.  The  first  of  all  the  stacks  to  detonate  is  called  the  initial 
donor.  BWACO  allows  the  user  to  specify  which  stacks  may  act  as  initial  donors.  It  is 
not  generally  possible  to  predict  which  ammunition  stack  will  act  as  an  initial  donor  when 
a  mass-detonation  accident  occurs.  Therefore,  in  most  applications,  the  user  should 
specify  that  all  stacks  may  act  as  initial  donors.  BWACO  will  issue  a  warning  whenever 
any  stacks  may  not  act  as  initial  donors.  (See  Input  Files.) 


Input  Files 


An  editor  that  produces  pure  ASCII  files  (without  hidden  characters),  such  as  the 
DOS  editor,  must  be  used  to  create  problem  input  files.  A  word  processing  program 
should  not  be  used  unless  the  file  can  be  saved  in  ASCII  format.  BWACO  looks  for  input 
in  a  file  with  the  extension  inp.  A  descriptive  file  name  (up  to  eight  characters)  may  be 
used.  The  file  must  conform  to  the  format  (see  Template)  shown  in  the  following 
example: 


Four  25-klb  Stacks,  Separation=200  m  (four.xnp) 


I — Xmin — I  I — Xmax — |  | 

0.  0. 

I - xs - 1  I - Ys - 1  I 

-100.  -100. 

100.  -100. 

-100.  100. 

100.  100. 


-Ymin —  |  |  — Ymax —  | 

0.  0. 

— Ws - II - ^Vs 1  D 

25000.  2500.  1 

25000.  2500.  1 

25000.  2500.  1 

25000.  2500.  1 
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The  first  line  of  the  file  contains  a  descriptive  problem  title  of  the  user’s  choice.  It  may 
contain  up  to  80  characters.  It  is  used  to  label  the  plots,  centered  in  the  plot  label  box 
at  the  top  of  both  the  Coalescence  and  Peak  Overpressure  Maps.  The  second  line  is  a 
ruler  line  delimiting  the  fields  in  which  the  bounds  of  the  region  of  interest  are  placed. 
These  bounds,  which  appear  on  the  third  line,  determine  the  region  in  which  the  stations 
are  distributed.  When  the  input  bounds  for  either  coordinate  axis  are  equal  (Xmin=Xmax 
or  Ymin=Ymax),  BWACO  computes  appropriate  bounds.  The  fourth  line  is  another  ruler 
line  specifying  the  fields  in  which  stack  data  is  placed.  Up  to  100  stacks  are  described, 
one  to  a  line,  on  the  remaining  lines.  The  stack  coordinates,  net  equivalent  explosive 
weight,  and  propagation  velocity  are  specified.  In  addition,  the  D  (initial  donor)  switch  is 
set  to  0  or  1  to  specify  whether  or  not  each  stack  may  act  as  an  initial  donor.  For  most 
problems,  the  initial  donor  switch  should  be  set  to  1  for  all  stacks.  (See  Initial  Donors.) 


Menus 

Menu  selections  are  made  by  entering  a  single-digit  number  at  the  prompt.  A 
menu  is  exited  by  entering  0  or  its  equivalent,  simply  striking  enter. 

The  Configuration  Menu  allows  selection  of  the  appropriate  video  display  and 
printer.  It  may  appear  in  any  of  the  following  three  forms; 

BWACO  CONFIGURATION  MENU 

Default  Display:  VGA  Mode  18  640x480/16 

Default  Printer:  Postscript  Printer 
■» 

Enter  1  to  change  the  default  display 
2  to  change  the  default  printer 
0  to  proceed  to  the  Main  Menu 


BWACO  CONFIGURATION  MENU 

Default  Display:  VGA  Mode  18  640x480/16 
Default  Printer:  Postscript  Printer 

Enter  1  to  change  the  default  display 
2  to  change  the  default  printer 
0  to  return  to  the  Main  Menu 


99 


BWACO  CONFIGURATION  MENU 

Default  Display:  VGA  Mode  18  640x480/16 
Default  Printer:  Postscript  Printer 

Enter  1  to  change  the  default  display 
2  to  change  the  default  printer 
0  to  return  to  the  Plotting  Menu 


The  Main  Menu  allows  selection  of  the  problem  execution,  plotting  or 
reconfiguration  functions,  as  available.  It  may  appear  in  any  of  the  following  three  forms: 

BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 
2  to  reconfigure  the  session 
0  to  end  the  session 


BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 

2  to  plot  the  current  problem 

3  to  reconfigure  the  session 
0  to  end  the  session 


BWACO  MAIN  MENU 

Enter  1  to  run  the  current  problem 

2  to  plot  the  previous  problem 

3  to  reconfigure  the  session 
0  to  end  the  session 


The  Plotting  Menu  allows  selection  of  Coalescence  or  Peak  Overpressure  Maps 
on  either  the  screen  or  the  printer.  It  appears  as  follows: 


BWACO  PLOTTING  MENU 


Enter  1 
2 

3 

4 

5 
0 


to  plot  the  Coalescence  Map  (Screen) 
to  plot  the  Coalescence  Map  (Printer) 
to  plot  the  Peak  Overpressure  Map  (Screen) 
to  plot  the  Peak  Overpressure  Map  (Printer) 
to  reconfigure  the  session 
to  return  to  the  Main  Menu 
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Peak  Overpressure  Map 

An  example  of  a  Peak  Overpressure  Map  is  shown  in  Figure  B-4.  Peak 
overpressure  contours  are  labeled  to  indicate  the  worst-case  values  at  each  point  within 
the  region  of  interest.  The  contour  levels,  which  are  listed  in  the  legend,  are  those 
specified  in  AR  385-64  for  the  protection  of  a  variety  of  assets.  For  comparison 
purposes,  an  additional  contour  representing  the  application  of  current  regulatory 
requirements  is  also  shown.  The  positions  and  equivalent  explosive  weights  of  the  stacks 
as  well  as  the  center  of  charge  and  total  equivalent  explosive  weight  are  also  shown  in 
the  legend.  (See  Coalescence  Map  and  Plotting  Distortion.) 


Plotting  Distortion 

Plotting  distortion  occurs  when  the  algorithm  used  to  extend  the  axes  to  a  natural 
termination  point  treats  the  X  and  Y  axes  differently.  To  eliminate  the  distortion,  it  is 
necessary  to  rerun  the  problem  after  modifying  the  input  file  to  specify  different  region 
limits.  (See  Coalescence  Map  and  Peak  Overpressure  Map.) 


Printers 

The  following  printers  are  supported: 

single-density  printer 
double-density  printer 
quad-density  printer 
Epson  printer  (plot  mode) 
Postscript  printer 
Hewlett-Packard  Laser  Jet 
Hewlett-Packard  Think  Jet 


Region  of  interest 

The  region  of  interest  is  the  region  covered  by  the  Coalescence  and  Peak 
Overpressure  Maps.  It  is  specified  by  giving  values  of  Xmin,  Xmax,  Ymin,  and  Ymax  in 
the  user’s  input  file.  If  the  user  specifies  values  of  Xmin  and  Xmax  which  are  equal, 
BWACO  will  compute  the  range  of  the  region  of  interest  in  the  X  direction.  An  analogous 
algorithm  applies  in  the  Y  direction.  (See  Stations.) 
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Stations 


When  running  a  problem,  BWACO  positions  3,600  stations  within  the  region  of 
interest.  These  are  not  generally  apparent  to  the  user.  Coalescence  is  checked  at  each 
of  these  stations  for  each  initial  donor.  (See  Region  of  Interest.) 


Template 

The  file  template .  inp  is  a  template  for  preparing  input  files.  It  can  be  copied 
as  the  first  step  in  preparing  a  new  input  file.  It  includes  the  two  ruler  lines  and  appears 
as  follows: 


I  — Xmin —  I  I  — Xmax —  |  |  — Ymin —  1  |  — Ymax —  | 

I - Xs - 1  I - Ys - 1  I - Ws - 1  I - Vs - 1  D 

(See  Input  Files.) 

Video  Displays 

The  supported  displays  are  listed  in  the  following  table: 


Type 

Mode 

Resolution 

Colors 

CGA 

4 

320x200 

4 

CGA 

5 

320x200 

4 

CGA 

6 

640x200 

BW 

Hercules 

720x348 

BW 

EGA 

13 

320x200 

16 

EGA 

14 

640x200 

16 

EGA 

15 

640x350 

BW 

EGA 

16 

640x350 

16 

VGA 

17 

640x480 

BW 

VGA 

18 

640x480 

16 

ATI 

81 

640x480 

16 

ATI 

82 

752x410 

16 

ATI 

83 

800x560 

16 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  ADMINISTRATOR  DTIC 
ATTN  DTIC  DDA 
CAMERON  STATION 
ALEXANDRIA  VA  22304-6145 

1  CDR  USAMC 
ATTN  AMCAM 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  DIR  USARL 

ATTN  AMSRL  OP  SD  TA 

RECORDS  MGMT 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIR  USARL 

ATTN  AMSRL  OP  SD  TL 
TECH  LIB 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DIR  USARL 

ATTN  AMSRL  OP  SD  TP 
TECH  PUBS  BR 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

2  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  TDC 
PICATINNY  ARSNL  NJ  07806-5000 

1  DIR  BENET  LABS 

ATTN  SMCAR  CCB  TL 
WATERVLIET  NY  12189^050 

1  DIR  USA  ADVANCED  SYSTEMS 

R&A  OFC 

ATTN  AMSAT  R  NR  MS  219  1 
AMES  RESEARCH  CENTER 
MOFFETT  FIELD  CA  94035-1000 

1  CDR  US  ARMY  MICOM 

ATTN  AMSMI  RD  CS  R  DOC 
REDSTONE  ARSNL  AL  35898-5010 

1  CDR  US  ARMY  TACOM 

ATTN  AMSTA  JSK 

ARMOR  ENG  BR 
WARREN  MI  48397-5000 


1  DIR 

USA  TRADOC  ANALYSIS  CMD 
ATTN  ATRC  WSR 
WSMR  NM  88002-5502 

1  CMDT 

US  ARMY  INFANTRY  SCHOOL 
ATTN  ATSH  CD  SECURITY  MGR 
FORT  BENNING  GA  31905-5660 


ABERDEEN  PROVING  GROUND 

2  DIR  USAMSAA 
ATTN  AMXSY  D 

AMXSY  MP  H  COHEN 

1  CDR  USATECOM 
ATTN  AMSTE  TC 

1  DIR  USAERDEC 
ATTN  SCBRD  RT 

1  CDR  USACBDCOM 
ATTN  AMSCB  CB 

1  DIR  USARL 

ATTN  AMSRL  SL  I 

5  DK  USARL 

ATTN  AMSRL  OP  AP  L 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  DIRECTOR 

US  ARMY  TECHNICAL  CENTER  FOR 
EXPLOSIVE  SAFETY 
ATTN:  SMCAC-ESL/C  DOYLE 
SAVANNA  IL  61074-9639 

3  US  ARMY  WATERWAYS  EXPERIMENT 

STATION 

ATTN:  CEWES-SE/ 

K  DAVIS 
MFORD 
C  JOACHIM 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  39180-6199 

3  CHAIRMAN 

DEPARTMENT  OF  DEFENSE  EXPLOSIVE 
SAFETY  BOARD 
ATTN:  DDESB-KT/ 

J  WARD  (2  CP) 

C  CANADA 

2461  EISENHOWER  AVE 
ALEXANDRIA  VA  22331-0600 

1  DEFENSE  AMMUNITION  LOGISTICS 

ACTIVITY 

ATTN:  AMCPM-AL/R  ROSSI 
PICATINNY  ARSENAL  NJ  07806-5000 

1  OFFICER  IN  CHARGE 
WHITE  OAK  LABORATORY 
NAVAL  SURFACE  WARFARE  CENTER 
ATTN:  R15/MSWISDAK 

10901  NEW  HAMSPHIRE  AVE 
SILVER  SPRING  MD  20902-5000 

2  NAVAL  FACILmES  ENGINEERING  SERVICES 

CENTER 

ATTN:  CODE  LSI/ 

J  TANCRETO 
RMURTHE 

PORT  HUENEME  CA  93043 

2  MAXWELL 

S-CUBED  DIVISION 
ATTN:  C  NEEDHAM 
S  HKIDA 

2501  YALE  BLVD  SE 
SUITE  300 

ALBUQUERQUE  NM  871064215 


3  SOUTHWEST  RESEARCH  INSTITUTE 

ATTN:  P  BOWLES 
C  OSWALD 
E  ESPARZA 
6220  CULEBRA  ROAD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

1  DIRECTOR 

EXPLOSIVES  HAZARD  REDUCTION 
DIRECTORATE 
ATTN:  JJENUS 
ASC/YOCO  (EHR) 

EGLDMAFBFL  32542 

2  US  ARMY  CORPS  OF  ENGINEERS 
ATTN:  CEHND-ED-CS/ 

R  WRIGHT 
RHASSE 
PO  BOX  1600 

HUNTSVILLE  AL  35807-4301 

3  COMMANDER 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  T  KENNEDY 
D  LINGER 
G  ULLRICH 

WASHINGTON  DC  20305-1000 

1  COMMANDER 

FIELD  COMMAND  DNA 
ATTN:  FCITS  (MCCRORY) 

KIRTLAND  AFB  NM  87117-5669 

1  HQDA  (SARD) 

ATTN:  MRMLARCHIK 
WASHINGTON  DC  20310-0001 

1  COMMANDER 

US  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
ATTN:  MONA-NU 
7150  HELLER  LOOP  SUITE  101 
SPRINGFIELD  VA  22150-3198 

1  ABERDEEN  RESEARCH  CENTER 

ATTN:  N  ETHRIDGE 

PO  BOX  548  . 

ABERDEEN  MD  21001 
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COPIES  ORGANIZATION 


1  CARPENTER  RESEARCH  CORPORATION 

ATTN:  J  CARPENTER 
PO  BOX  2490 

ROLLING  HEUS  ESTATES  CA  90274 

1  LOGICON  R&D  ASSOCIATES 

ATTN:  GGANONG 
PO  BOX  9335 

ALBUQUERQUE  NM  87119 

1  TECH  REPS  INC 

ATTN:  FMCMULLAN 
5000  MARBLE  NE  SUITE  222 
ALBUQUERQUE  NM  87110 

ABERDEEN  PROVING  GROUND 

20  DIR,  USARL 

ATTN:  AMSRL-WT-T/W  MORRISON 
AMSRL-WT-TB/ 

RFREY 
O  LYMAN 
T  WATSON 
F  GREGORY 
V  BOYLE 
W  HILLSTROM 
E  MCDOUGAL 
W  LAWRENCE 
,  K  BENJAMIN 
T  DORSEY 
AMSRL-WT-TD/ 

A  DAS  GUPTA 
J  SANTIAGO 
AMSRL-WT-NC/ 

R  LOTERRO 
R  RALEY 
J  SULLIVAN 
S  SCHRAML 
AMSRL-CI-CA/ 

A  CELMDSrS 
C  ZOLTANI 
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1  SNPE 

ATTN:  LALLAIN 

BP2  91710  VERT-LE-PETIT 

FRANCE 

2  MINISTRY  OF  DEFENSE 
EXPLOSIVES  STORAGE  AND 

TRANSPORTATION  COMMITTEE 
ATTN:  J  CONNOR 
M GOULD 

ROOM  1818/EMPRESS  STATE  BUILDING 

LONDON  SW6  ITR 

ENGLAND 

1  AGENCY  FOR  DEFENSE  DEVELOPMENT 

ATTN;  S  -  Y  SONG 
PO  BOX  35 

YUSEONG  TAEJON  305-600 
REPUBLIC  OF  KOREA 

1  AUSTRALIAN  ORDNANCE  COUNCIL  DOD 

ATTN;  JGOOLD 
CANBERRA  ACT  2600 
AUSTRALIA 

1  NATIONAL  DEFENSE  RESEARCH 

ESTABLISHMENT  FOA 
ATTN:  RFORSEN 
PO  BOX  551 
TUMBA  S- 14725 
SWEDEN 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-645 _ Date  of  Report  December  1994 _ 

2.  Date  Report  Received _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) . . . 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ _ _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 
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